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ABSTRACT 


The understanding of liquid fuel droplet combustion is very important for 
designing and developing low emission and energy efficient combustion systems. 
The present thesis is devoted to investigate the combustion aspects of a single 
isolated droplet nmnerically and experimentally. Droplet combustion experiments 
for both single and multi component fuels at normal pressure and temperature are 
performed under normal gravity condition. The cf law diameter regression rate 
behavior is observed for both single and multicomponent fuel. The burning rate 
constant for multicomponent fuel droplet combustion lies between that for its 
individual constituents. The droplet combustion is modeled by solving one 
dimensional, unsteady governing equations numerically. In this CFD code, 
variable properties, gray gas model for radiation heat transfer are included. The 
effects of ambient oxygen level on droplet combustion are very important to study 
as it varies in the course of combustion. Burning rate constant as well as the flame 
diameter increases with the increase in ambient oxygen level. The flame 
temperature is increased leading to gasification of fuel droplet. The inclusion of 
gray gas model for radiation heat transfer improves the prediction of droplet 
burning rate and flame temperature. It was found out that radiation heat transfer 
(heat loss at low ambient temperature and heat gain at high ambient temperature) 
for bigger size droplet is more predominant compared to smaller droplet. As a 
result opposite influence of the initial droplet diameter on the droplet combustion 
can be observed between the higher and lower temperature. The burning rate 
constant with increase in initial droplet diameters is increased and decreased at 
high and low ambient temperature respectively. In the third phase, the 
multicomponent fuel droplet combustion is investigated by using distillation curve 
model. The transient analysis of the surface properties with different fuel 
composition has also been carried out extensively in the study. The data presented 
in the present study matches well with the available previous study. It was found 
out that multicomponent fuel droplet combustion has three dominant periods, 2 
law periods separated by transition period with a gradual change. 
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Chapter 1 

INTRODUCTION 


1.1. Background 

The first step in human civilization is the invention of fire. The further mastery over 
combustion and its applications played most vital role to drive the human 
civilizations. Industrial era was started with the invention of steam engine by James 
Watt. This was the first step to convert the heat of combustion to the mechanical 
power. Even today 90% of the energy used by the human being is being obtained by 
burning of fuels such as, coal, petroleum oils, biomass, etc. The phenomenon of 
combustion finds applications in several modem systems such as internal 
combustion engine used in most of the automobiles, gas turbine engines in aircraft 
and power plant, boiler used for power generation and lots of other applications. 

However, the over exploitation of combustion process to meet the energy 
demands of modem man has created unmitigated problems of emissions which is a 
threat to very existence of life on the earth. The air pollutants such as carbon 
monoxide, NOx, SOx, unbumt hydrocarbons, etc. are generated due to incomplete 
combustion of fuel. Combustion process also produces large amount of smoke, ash, 
odors, noxious and benign gases. These gases have damage the environment around 
to a large extent for which there is a global concern to reduce the emission from the 
combustion systems. Beside this there is scarcity of fossil fuel across the globe. 
Hence, it is important to utilize the available fossil fuel and enhance the efficiency 
of the combustion systems. For which, there is an urgent need to understand the 
processes involved in the combustion and develop energy efficient and low 
emission combustion systems. 

Generally the liquid fuel is used in most of the combustion devices due to 
easy storage and high energy density. However, it can not be burnt itself efficiently 
like any other gaseous fuel. In order to increase the burning rate the contact area 
between the liquid surface and the ambient oxygen has to be increased, which can 
be possible by disintegrating the liquid fuel into smaller droplets. Disintegrating a 
bigger droplet into n numbers of smaller droplet total surface area can be increased 
by n times. Hence the ftiel can be burnt time faster than the previous case. For 
example, burning time can be enhanced by 10‘ times to disintegrate a 3 mm 
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diameter bigger droplet into 10^ numbers of 300 pm diameter droplets. In any liquid 
fuel combustion system, fuel is atomized into spray containing large numbers of 
droplets with different diameters. But the understanding of spray combustion is 
impeded by the complex processes involved in it. However the whole spray 
consisting of a wide range of droplet diameters can be well represented by a single 
droplet diameter generally in terms of Sutter Mean Diameter (SMD). The SMD of a 
spray is used as an important parameter during the design of a combustor. Hence 
the study of the droplet combustion is very important to understand the physical 
mechanism during its combustion process. 

1.2. Droplet Combustion 



Fig. 1.1: Liquid fuel droplet combustion in (a) Micro gravity (b) Normal gravity (c) 
Forced convection conditions [47] 

The processes involved in the combustion of a single droplet rmder three different 
conditions such as (a) Microgravity, (b) Normal Gravity and (c) Forced Flow are 
described in Fig. 1.1. During droplet combustion, liquid fuel vaporizes and mixes 
vsdth the air due to mass diffusion. As the temperature outside the droplet is higher 
than the self ignition temperature of the fuel, combustion takes place. Heat and 
products are generated, which are transferred due to heat and mass diffusion 
respectively. When the droplet is burnt under normal gravity condition, the high 
temperature product gas moves upward due to buoyancy force and natural 
convection heat transfer takes place, as shown in Fig. 1.1(b). The burning rate is 
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enhanced considerably due to convection heat transfer. However under zero or 
micro gravity conditions the flame surface becomes almost spherosymmetric as 
shown in fig. 1.1(a). For this case the ratio of the flame radius to the droplet radius 
is defined as flame stand-off ratio. When the flow takes place over the droplet 
surface as shown in Fig. 1.1(c), the flame is aligned along the free stream velocity. 

Classical Theory on Droplet Combustion: 

For simplicity, droplet burning is considered as follows in the classical theory of 
droplet combustion; 

(1) Spherical symmetry; (2) No spray effect; (3) Combustion rate controlled 
by mass transfer to the flame due to diffusion; (4) Constant pressure process; (5) 
Thin flame approximation; (6) Constant gas phase properties; (7) Quasi steady 
process assumption; (8) Constant and uniform droplet temperature; (9) Neglecting 
Soret, Dufour and radiation effect, (10) unity Lewis number, (1 1) No radiation heat 
loss. 

By using this approximation the instantaneous droplet diameter can be 
represented by the initial droplet diameter and the burning time as follows; 

d^=dl-K,t ( 2 . 1 ) 

where, d and d^ are the instantaneous droplet diameter at time t and initial droplet 
diameter. Kf, is the burning constant (also known as effective burning rate in 
literature) which can be expressed as 

K^=?±J-£^\n(\ + B) ( 2 . 2 ) 

dt p, 

where, and -are the gas and liquid phase density, is the gas phase thermal 
diffusivity at the surface and B is the Spalding transfer number. And the total 

burning time can be represented zs . 
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1 .3. Motivation for the Present Work 

In any combustion device, the performance and size of combustor are strongly 
dependent on the total allowable time provided to bum the droplet. If the droplet 
does not get the required time to bum, unbumt hydrocarbons are exhausted out to 
the atmosphere. As a result the efficiency of the combustion systems is reduced 
drastically while releasing large amount of pollutants to the atmosphere. Most of the 
time, the liquid fuel droplet are burnt under high pressure and temperature. Beside 
this the concentration of oxygen level in the combustor may get changed with time, 
which affects its burning efficiency. Hence it is very important to study the effects 
of temperature, effects of oxygen concentration on the droplet burning, which is not 
being addressed adequately in the literature. 

Most of the practical liquid fossil fuel generally consists of several blends of 
hydrocarbons with different physical and chemical properties. However, the 
literature review reveals that the understanding of the multicomponent fuel droplet 
combustion is at infancy stage. Hence, in this thesis an attempt is made to study the 
multicomponent fuel droplet combustion by using both experimental and numerical 
tools. 

1.4. Objective 

The main objectives of the current work are the followings. 

i) Design and develop an experimental setup to investigate the single isolated 
droplet combustion under normal gravity. 

ii) To study the single and multi component fuel droplet combustion under 
normal atmospheric pressure. 

iii) To develop unsteady, 1-D, variable properties, single step chemistry, 
numerical model for droplet combustion under microgravity for single and 
multi component fuel droplet combustion. 

iv) To study the transient variation of droplet diameter square, flame structure, 
surface parameters such as surface temperature gasification rate, etc. 

v) Effect of ambient oxygen concentration on the single component fuel 

droplet combustion. 

vi) Study the effects of radiation heat transfer on the droplet combustion. 

vii) Study the inverse influences of the initial diameter in high and low ambient 
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temperature. 

viii) Study of multicomponent fuel droplet combustion. Changing in surface 
parameters during multicomponent fuel droplet combustion. 

1 .5. Organization of the Thesis 

The background, motivation and objective of this research are contained in the 
present chapter. Previous research on the field of droplet combustion is discussed in 
brief in chapter 2. The major previous experimental and numerical works are 
reviewed critically in this chapter reporting the anomalies present in the literature. 
Based on the literature survey the objectives of the present thesis are set out. The 
experimental setup and the results of the present work are discussed in chapter 3. In 
chapter 4, the mathematical model for this analysis is elaborated, which is followed 
by the numerical results and discussion in chapter 5. Results and discussions are 
presented in two sections Section 5.1 containing the simulation results for single 
component fuel droplet combustion followed by the results for the multicomponent 
fuel droplet combustion in Section 5.2. The effects of initial oxygen concentrations, 
effects of radiations are discussed extensively for single component fuel droplet 
combustion. However, the effects of compositions on the transient variation of 
different burning parameters are being discussed for multicomponent fuel droplet 
combustion. The last chapter (Chapter 6) contains the conclusion and the 
recommendations for the further studies. 

1.6. Summary 

In the present chapter a brief insight into the basics of the droplet combustion is 
presented. It is quite clear from the study that droplet combustion is very important 
to study. The processes involved during the combustion of single isolated droplet 
under three different conditions such as, normal gravity, micro-gravity and forced 
flow conditions are being discussed briefly. The classical theory of droplet 
combustion is enumerated briefly highlighting the universality of dr law. The 
motivation of the present work is described mentioning the current issues to be 
addressed in this thesis. The objectives of the present thesis are enumerated 

adequately. 
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Chapter 2 

LITERATURE REVIEW 

2.1. Classical Theories on Liquid Fuel Droplet Combustion 

First successful and remarkable investigations of droplet combustion were carried 
out theoretically by G A E Godsave [18] and D B Spalding [68] independently in 
1953. Godsave [18] presented two distinct models such as (i) Low Temperature 
Evaporation and (ii) High Temperature Evaporation and Combustion. He found out 
that Low Temperature Evaporation is more likely to occur when the temperature 
difference between the droplet and the ambient is less. On the other hand, the High 
1 emperature Evaporation and Combustion can occur when the difference in 
temperature between the droplet and surround atmosphere is considerable. He 
concluded that the volatility of the fuel, expressed in terms of vapor pressure, does 
not have a predominant effect on the burning rate. However the instantaneous mass 
rate of burning is proportional to its diameter. 

Spalding [68] was first to introduce a non-dimensional transfer number, B in 
his quasi steady state model that controls the mass burning rate. Subsequently this 
was leveled as Spalding’s transfer number. The mass transfer from liquid fuel 
droplet surface depends only on Spalding’s transfer no, B. In this model, the rate 
controlling phenomenons were found to be dependent on the physical properties, 
heat and mass transfer. However the chemical kinetics associated with combustion 
does not play any critical role. 

2.2. Experimental Studies 

Experimental studies on droplet combustion was pioneered by Kumagai, S., and 
Isoda, H. [28] in 1955. Subsequently they investigated the combustion of a single 
isolated droplet in microgravity successfully in 1957 [30]. In that experiment they 
studied the effects of gravity on the droplet combustion in terms of flame structure 
and the droplet burning rate. ■ 

In 1971, Kumagai, Sakai and Okajima [31] performed the first unsupported 
droplet combustion experiments using a 0.8 s drop tower. In these experiment, the 
burning rate of n-heptane was reported to be nearly constant, 0.78 mm /s. 
Considerable numbers of experiments supported their findings (Yang and Avedisian 
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[80]; Hara and Kumagai [20]). Hence the burning rate constant for an isolated n- 
heptane droplet was remained as a fundamental parameter for more than two 
decades. Later on, advanced theories supported by experimental data revealed the 
earlier conclusion not to be right. 

2.3. Parametric Studies of Droplet Combustion 

The droplet combustion at high pressure was modeled by Rosner [60], 
Chervinsky [7]. They found out that burning rate is proportional to Faeth et al. 
[13] also conducted supercritical combustion studies up to pressure of 100 atm., 
under micro-gravity conditions. Their experimental results supported the previous 
modeling results. Later on the transient variation of droplet combustion by varying 
operating parameters for a wide range were investigated extensively by several 
research groups such as, Hubbard et al. [24], Law [36], Law and Sirignano [37], 
Saitoh and Nagano [61]. Among all the results of Hubbard et al. [24] are very closer 
to the existing experimental data. Prakash and Sirignano [57], [58] modeled the 
internal circulation in the droplet due to shear generated at its surface depicting the 
liquid and air interaction. 

Egawa et al. [12] examined the pressure effects on combustion of fuel 
droplet under zero gravity. They concluded that the evaporation constant increases 
with increase in pressure and its effects on the evaporation constant are very small. 
Flame stand-off ratio decreases with pressure and its variation carmot be explained 
by the classical quasi-steady theory. Subsequently Sato et al. [64] carried out 
experiments on the effects of natural convection on high pressure droplet 
combustion. They studied the burning of n-octane fuel droplet under both normal 
and microgravity conditions. They concluded that the burning rate increases in the 
sub-critical zone and decreased in the super-critical zone in both micro and normal 
gravity condition with the increase in pressure and attains a maximum value at the 
critical pressure. The burning rate increases with the natural convection and its 
effects are relatively pronounced at high ambient pressure. 

The effect of acoustical excitation on the droplet evaporation and 
combustion was studied experimentally by Saito et al. [63] in 1995. They concluded 
that the acoustical excitation promotes the mixing between the fuel vapor and the 
surrounding gases leading to augmentation in heat and mass transfer during the 
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combustion process. As a result the burning rate was enhanced. Ghenai et al. [15] 
extended it further and concluded that burning rate increases more that 20% when 
the droplet is situated at the pressure node (velocity anti-node) of the acoustic wave. 
In contrast no such effects were observed at velocity node (pressure anti-node). 

2.4. Effects of Radiation on Droplet Burning 

Radiations heat losses from the flame surface are two types: luminous radiation and 
non-luminous radiation. Non-luminous radiations are mainly associated with the 
combustion gases, such as CO 2 and H 2 O. In contrast, oxygen and nitrogen are 
considered as transparent to the thermal radiations. Carbon particles or soot are 
responsible for the luminous radiation (Tien and Lee [70]). 

'I’he classical quasi-steady theory does not take account of the radiation 
effects on the droplet combustion. However, the effect of radiation was analyzed by 
Godsave [18] by considering the heat directly transferred from hot flame zone to the 
cold droplet surface. The reflection and radiative losses from the droplet was not 
considered in this model. Total heat gain due to radiation at droplet was 
approximated as, = Aa^a ; where, A is the drop surface area, is the 
absorptivity of drop surface, cty^is the absorptivity of flame front and ctj^ is the 

Stefan-Boltzman’s constant. In this analysis, the flame temperature, mean droplet 
surface were assumed to be 2000 K, 323 K. Beside this they also approximated the 
flame surface and drop surface absorptivity to be 0.03 and 0.05 respectively during 
the burning of 1 .42 mm diameter benzene droplet. In his study it was concluded that 
20% of heat is transferred at droplet surface due to radiation and absorbed radiation 
is proportional to the diameter of the droplet. It was suggested that the radiation 
effect can be neglected for smaller droplets. 

In another study. Beer and Howarth [5] concluded that any small gas 
volume in any thermodynamic system must be in a thermodynamic equilibrium. 
The product species such as H 2 O and CO 2 are the considered to be the major 
contributors for radiative heat transfer. According to the study of Abdel-Khalik et 
al. [1], 40% of total heat is transferred due to radiation. Out of which 30% is 
contributed by the soot radiation, while 70% is due to gas phase radiation. In their 
study H 2 O and CO 2 were considered as only the participating gases. 

Currently, two types of radiation models such as Gray gas and Non-gray gas 
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model are being used in the literatures. In Gray gas model the radiation properties 
are assumed to be remained constant with the wavelengths. However, the variations 
of radiation properties with the wave-lengths are considered in non-gray gas narrow 
and wide band models. In the weighted sum gray gas model, the gas is represented 
by a number of gray gases with different absorption coefficients. 

In most of the studies such as, Saitoh et al. [62], Marchese and Dryer [43], 
Marchese et al. [44], Kumar et al. [33], Kazakov et al. [27], gray gas modeling 
originally pioneered by Viskanta and Merriam [75] had been used. Saitoh et al. [62] 
used a constant Planck mean absorption coefficient equal to 0.35 m'‘. While 
Marchese and Dryer [43], Marchese et al. [44], Kazakov et al. [27], Kumar et al. 
[33] had employed the variable Planck mean absorption coefficients, which was 
evaluated locally based on the partial pressure weighted sum of mean absorption 
coefficient for each participating gas at the local grid temperature. 

In 1992 a comprehensive mathematical model incorporating multi- 
component molecular diffusion, complex chemistry, and heterogeneous process to 
investigate a variety of chemically reacting flow problem was introduced by Cho et 
al. [10]. This model emphasized the elementary chemical and physical process 
rather than complex flow geometry. They concluded that extinction of smaller 
droplets is caused by purely finite-kinetics, while extinction of larger droplets is 
caused by excessive radiative heat loss which eventually limited the kinetics. 

Effects of non-luminous thermal radiation in micro-gravity droplet 
combustion were examined by Marchese and Dryer [43]. They numerically 
simulated the transient combustion of methanol-water mixture droplets. Their 
results also indicate that the effect of radiative heat loss cannot be neglected for 
larger initial droplet diameters. 

Transient combustion of spherical single and multicomponent n-Alkane 
droplets with non-luminous gas phase radiative heat transfer with complex chemical 
kinetics was numerically simulated by Marchese et al. [44]. By including radiation 
in this model, a significant decrease in maximum flame temperature was observed. 
The normalized flame position decreases as a consequence of radiative losses. With 
increase in initial diameter, their model predicts lower burning rate and eventually 
leading to radiative extinction. 

Baek et al [4] investigated the influence of non-gray radiation on the droplet 
combustion using weighted sum of gray gas model assuming CO 2 and H 2 O to be 
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the major contributors. By using a constant Planck mean absorption coefficient of 
2.0 m they found out a good agreement between the results of gray-gas model and 
non gray gas model. 

lo investigate the effect of initial droplet diameter on combustion at 
different ambient temperature, Xu et al. [79] performed their experimental analysis 
of a single isolated n-heptane droplet at 298 K, 633 K and 1093 K. Their result 
predicted the inverse influences of initial droplet diameter in cold and hot 
ambiences. The burning rate respectively decrease and increase with the increase in 
droplet diameter for former and later case. The differences between the radiation 
heat gain by the droplet and heat loss to the ambience was suggested to be 
responsible for this. Subsequently this phenomenon at room temperature was 
reported in the numerical and experimental analysis of Kazakov et al. [27]. 

2.5. Droplets Interaction 

In actual combustion systems, droplet never bums in an isolated manner without 
any interaction among the adjacent droplets. Under this situation, the droplet- 
droplet (and flame-flame) interaction will affect the flame shape, flame dynamics, 
burning rates and extinction behaviors. Early investigations using linear arrays of 
droplets were performed by Miyasaka and Law [48], Xiong et al. [78]. Their results 
indicated that the average burning rate for the droplet arrays is lower than that of an 
isolated droplet. This phenomenon is attributed to oxygen starvation. Further 
studies indicated that burning rate increases with the decrease in separation 
parameter, which is defined by l/do, where I is the distance between two droplet 
centre and do is the initial droplet diameter. Recently, the study performed by 
Mikami et al. [46] concluded that for both micro gravity and normal gravity, there 
are reductions in the burning time (increase in burning rate) as the separation 
parameter; l/do is decreased from infinity (corresponding to isolated droplet 
combustion). For normal gravity results, as the separation distance is reduced, 
oxygen transport to the flame is enhanced due to natural convection around each 
droplet. However additional reduction in separation distance is resulted in oxygen 
starvation, which increases the burning time (decrease in burning rate) drastically. 
In the case of microgravity, the radiation heat gain at the droplet surface initially 
increases when l/do is decreased from infinity resulting in the decrease in burning 
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time. Further decrease in l/do creates the oxygen starvation effect and burning time 
increase dramatically. Explanation of such behavior of droplet burning rate for 
micro-gravity was provided by Okai et al. [51] in the experimental result of droplet 
near wall. 

2.6. Multicomponent Fuel Droplet Combustion 

Most of the fuels used in practical combustors are multi-component in nature. The 
processes involved in the burning of multicomponent fuel droplet are much more 
complex than single component fuel combustion which is mainly due to liquid 
phase mass transfer. The phase change characteristics of a multicomponent fuel 
depend on the extent of miscibility among its constituents. Several models for the 
multi-component fuel droplet combustion are being put forwarded by researchers. 
Distillation curve model, diffusion limit model and shell model sure some of the 
important models for liquid phase mass transfer during multicomponent fuel droplet 
combustion. In diffusion limit model, liquid phase transport is considered to be due 
to the mass diffusion. It is generally assumed that there is no any internal circulation 
present inside the liquid phase. In contrast, liquid phase transport is considered to 
occur due to internal circulation in distillation curve model resulting in uniform 
distribution of fuel constituents inside the droplet. Hence, the species concentrations 
are assumed not to be varied spatially. The rapid mixing or distillation curve model 
of multi-component droplet burning was initially investigated by Faeth [14], and 
Law [36]. In 1976 Law [36] introduced the shell model assuming that single- 
component fuel vaporizes in each shell. By further simplifications he obtained the 
results which were closely approximated to the bulk vaporization behavior as 
described in the distillation curve model. 

Law and Law [40] developed a model for asymptotic analysis on diffusion 
limited quasi-steady regime assuming the extremely slow rate of liquid phase mass 
diffusion. No changes in the temperature and mass fraction profile on the liquid 
phase were considered in their study. Their results predicted that the fraction of 
mass flux leaving the droplet surface for any species was equal to the initial liquid 
mass fraction of that species and (f law prevails during the burning of 
multicomponent fuel droplet combustion. Experimental results obtained for 
alchohal-alkane binary droplet combustion by them confirm the law behavior. 



Results presented for 50% Undecane and 50% Octanol mixture is shown in Fig.2.1. 



t/d^ t (sec) 

0 


Fig. 2.1: Experimental results for (a) the existence of d^ for multicomponent fuel 
droplet combustion by Law and Law [40], (b) temporal variation of the droplet 
diameter square with time for heptane/hexadecane droplet by Wang et al. [76]. 

Combustion and Micro-explosion studies on a freely falling multicomponent 
fuel droplets carried out by Wang et al. [76]. The combustion characteristics of 
isolated, low Reynolds number, multicomponent droplets falling in a hot, oxidizing 
gas flow were studied. Three staged combustion behavior were reported and the 
mass diffusion were reported as dominant liquid phase transport mechanism. It was 
concluded that, during multicomponent fuel droplet combustion there exists an 
initial period when the volatile components in the surface layer gets vaporized 
droplet temperature remain relatively cold. Followed by a transition period when 
droplet surface temperature increases rapidly. In this period the burning rate is 
extremely low, the size shrinks and the temperature diminishes. Finally a quasi- 
steady period sustains during which the droplet interior concentration distribution 
remains almost constant and the surface region is more concentrated with the less 
volatile component. The results of the heptane/hexadecane droplet combustion 
presented in their study are shown in Fig. 2.1 (b). The results for the micro- 
explosion of a multicomponent fuel droplet consisting 60% propanol, 30% 
hexadecane, and 10% surfactant during the transitional droplet heating period is 
shown in Fig. 2.2 (b). It may happen during burning of droplet when less volatile 
constituents are at outside and the more volatile constituents are at the center shells. 
Thus the more volatile components superheat and try to come outside fast due to 
extra pressure generated during superheat. This may start homogeneous nucleation, 
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whose extremely fast gasification rate may lead to sufficient pressure build up and 
thus the droplet may eventually fragment. 
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Fig. 2.2: Experimental measurement of the square of droplet diameter vs. time for 

(a) hexadecane/heptane droplet burning in 1 atm. air (Yang and Avedisian, [80]); 

(b) A droplet consisting of 60% propanol, 30% hexadecane, and 10% surfactant 
(C.H. Wang, X.Q. Liu and C.K.Law, [76]). 


Tong and Sirignano [71] modeled the multicomponent transient droplet 
vaporization with internal circulation. The liquid phase internal circulation, 
transient droplet heating and species diffusion were taken care in their model. 

In the study of the major controlling parameters in the gasification behavior 
of multicomponent droplets Makino and Law [41] in 1988 concluded that the liquid 
phase Lewis number is not the appropriate parameter to indicate the influence of 
liquid phase diffusional resistance on the droplet gasification behavior. However, it 
was found out that Peclet number, which is defined as the ratio of a reference 
surface regression rate to the liquid phase mass diffusivity, have a crucial role in 
controlling the gasification from the liquid fuel surface. 

The experiment of unsupported hexane and hexadecane mixture under micro 
gravity was performed by Yang and Avedisian [80] using 1.2 s drop tower. Their 
results are plotted in Fig. 2.2 (a). Their data concluded that the evaluation of droplet 
diameter was similar to that of the classical cf law for a single droplet. In order to 
investigate the law behavior of multicomponent fuel droplet conclusion. Law and 
Law [40] concluded the similar nature of cf law. 
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To resolve the above anomalies between the results, predicting the nature of 
the burning history (cf law) for multicomponent fuel droplet, Shaw [66] extended 
this study of Law and Law [40]. His study on the multicomponent fuel droplet 
concluded that the droplet species profile change is slow relative to the change of 
droplet size. Hence, (f law surface regression closely followed. The mass flux 
fraction of a given species of the droplet surface is approximately equal to the 
volume averaged mass fraction of that species within the droplet. The ratio of the 
liquid phase species diffusion coefficient can be treated as a small parameter in 
contrast to the birming rate constant. 

The evolution of liquid-phase mass fraction and temperature profiles during 
combustion of fuel droplets composed of binary miscible mixtures of low volatility 
and high volatility constituents was studied by Shaw and Williams [67]. They 
assumed the small initial mass fractions of low volatility material or slow 
introduction of low volatility material into the liquid phase by absorption from the 
gas phase. The gas phase was considered to remain quasi-steady and the liquid 
temperatures were considered to be spatially uniform. Their result shows that at the 
droplet surface a boundary layer was developed in which the mass fraction of the 
low volatility components were increases with time. Both the law and non-^f^ law 
behaviors of the droplet surface regression were discussed. Their study concluded 
that for small initial values of the initial impurity mass fraction, the surface 
regression approximately follows the d‘ law nature until the mass fraction of the 
impurity reaches nem unity. Meanwhile, for sufficiently large values of the initial 
impurity, their numerical results demonstrated appreciable departures from the 
law behavior. 

A fully transient numerical study for multicomponent liquid fuel droplet was 
carried out by Mawid and Agarwal [45]. This numerical result demonstrate that due 
to the process of fuel vapor accumulation within the inner region to the flame and to 
the high liquid-phase mass diffusional resistance, a situation is reached where the 
concentration of fuel vapor in the inner zone of the flame or in the immediate 
vicinity of the droplet surface exceeds that at the droplet surface for certain values 
of the liquid phase Lewis number.' This result gives negative vaporization rate and 
no indication of flame extinction. Typical variation of fiP with time for 50-50 
heptane-decane mixture as reported by them is shown in fig. 2.3 (a). 
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Fig. 2.3: (a) Temporal variation of the droplet radius squared for the base case. The 
asterisk indicates the time when suppression of vaporization occurs (Mawid and 
Aggarwal [45]). (b) Predicted square of the droplet diameter vs. time for law, 
distillation curve and diffusion limited model (Glassman [17]). 

The results predicted by Glassman [17] are shown in Fig. 2.3 (b). The heat 
up period predicted by the rapid mixing model is longer since more heat is required 
to raise the temperature of the entire droplet (relative to the diffusion limit 
model).To study the effect of water concentration on flame luminosity and sooting 
tendency and disruptive burning, Jackson and Avedisian, [25] performed the 
experimental study of combustion of unsupported n-heptane droplets emulsified 
with water in microgravity. They concluded that the sooting is minimized for 
emulsions relative to pure heptane. A period was reported near the end of burning 
during which the droplet diameter is nearly constant, followed by disruptive 
burning events. 

Mandal [42] studied the combustion behavior of single isolated multi 
component fuel droplet using finite element method. He used adaptive grid 
generation method and variable gas and liquid phase properties. A highly efficient 
semi-analytical model was also developed. Both the models exhibit three dominant 
period of droplet burning, two d^-law periods separated by a transition period which 
exhibits a gradual change. In both the models, it was considered that the liquid 
phase mass transport occur only due to diffusion process and internal circulation, 
which is more likely to be occurred in the actual situation, was not considered.. 

Very recently in 2003, droplet evaporation study by using distillation curve 
model for kerosene fuel at normal and elevated pressure was performed by Burger 
et al. [6]. Their results were in good agreement to the diffusion limit model at 
different pressure. 
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2.6.1. Comparison of Results Available in Literature 

In the earlier section the investigations on the multicomponent fuel droplet 
combustion available in the open literature were discussed extensively. This works 
can be categorized broadly into two groups such as (i) the experimental and 
numerical studies presented by Law and Law [40], Yang and Avedisian [80], and 
Shaw [66], those predicts the classical ^ law and (ii) Wang et al. [76], Mawid and 
Aggarwal [45], those have reported the staged behavior of (f law. On the critical 
examination of the experimental data available in the literature, it was found out 
that there is large discrepancies exist in the data by various researchers even for the 
almost same composition of fuels. The variation of droplet diameter square with 
time reported by Wang et al. [76] and Yang and Avedisian [80] for heptane 
hexadecane mixture are compiled and plotted in Fig. 2.4 (a). It can be observed that 
the results presented by Wang et al. [76] exhibited non law behavior. In contrast 
the results of Yang & Avedisian [80] departed from the data of Wang et al. [76] and 
while retaining the classical (f law behavior. 



time (sec) t/d^ (sec/mm^) 

Fig. 2.4: Comparison of results: (a) The results predicted by Yang & Avedisian 
[80] and Wang et al. [76] for burning of heptane-hexadecane fuel mixture, (b) 
Results presented by Wang et al. [76]. 

It can be noted that Wang et al. [76] claimed the experimental data to follow 
staged (f law behavior. However it is conjectured that this changed in slope of 
diameter square with time may be due to experimental error. In order to asses their ^ 
results critically all the data are compiled and presented in the standard normalized 
form {d^/do vs. 't/do) in Fig. 2.4 (b). It can be observed that their data for 
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multicomponent fuel droplet combustion for two different compositions exhibit 
both the classical (f law and non law behavior. The diameter square change with 
time for multicomponent fuel droplet with composition of 20 % heptane and 80 % 
hexadecane follows the classical law. In contrast multicomponent fuel droplet 
with composition of 70 % heptane and 30 % hexadecane departs slightly from the 
classical law. However the results presented by Yang & Avedisian [80] for 
almost same composition of fuels depict cf law behavior. 

2.7. Conclusions 

The most important step in the study of droplet combustion was the development of 
the classical quasi-steady theory. Initial droplet combustion research efforts were 
directed towards ascertaining the validity of the classical quasi-steady theory. 
Though the (f law predicted by the classical theory can corroborate to the 
experimental data, but the flame stand off distance predicted by the classical quasi- 
steady theory departs substantially from that of the experimental data. The addition 
of radiation heat transfer effects in the theoretical model predicts the flame stand off 
ratio reasonably. Although several radiation heat transfer models have been 
developed the droplet burning rate, but several discrepancies exists in the literatures 
which needs to be addressed. Recently the experimental data reported by Xu et al. 
[79] revealed an interesting feature of opposite influence of the initial diameter to 
the burning rate constant particularly at high temperature. And this behavior of 
droplet burning has not been investigated numerically to best of our knowledge. 
Beside this, the effect of eimbient oxygen concentration on the single component 
fuel droplet combustion has not been addressed adequately, which will be 
investigated in detail in the present thesis. 

Several numerical and experimental studies on the multicomponent fuel 
droplet combustion have been reviewed in the present chapter. It was found out that 
three major models such as distillation curve model, diffusion limit model and shell 
model for liquid phase mass transport have been put forwarded by the researchers. 
It seems that the results of distillation curve model in comparison to the other model 
are closer to the experimental data. Hence, there is a need to develop a better model 
for multicomponent fuel droplet combustion. The surface parameters during 
multicomponent fuel droplet combustion are very important to study. 
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Chapter 3 

EXPERIMENTAL STUDIES 

3.1. Introduction 

In order to study the droplet combustion, an experimental setup has been deigned 
and developed from the scratch in the lab. It can withstand maximum pressure of 5 
MPa. The design calculation was done with the help of ASME Boiler and Pressure 
Vessel Design Code [3]. In this case the factor of safety of 3 is employed. Stainless 
steel is chosen as material for combustion chamber. The detailed description of the 
experimental setup and the experimental results are discussed below. 

3.2. Experimental Setup: 


Pressure Gauge Valve 



Fuel Injector 



Fig. 3.1: Schematic of the experimental setup (a) Front view (b) Side view. 


Two orthogonal views (Front and Side view) of experimental setup used for the 
present study are shown schematically in Fig. 3.1. The stainless steel box of 
dimension, 160 mmx 160 mmx 180 mm is used as a combustion chamber. All the 
walls of the combustion chamber are welded except the top wall, which is fixed by 
nut and bolt. The three side walls of the combustion chamber are equipped with 
circular glass windows in such a way that the droplet remains along the optical axis 
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of these windows during experimentation. A glass fiber of diameter 0.16 mm is 
used to suspend the droplet. This fiber is attached to the top wall with a specially 
designed holder. In order to hold a droplet at the tip of the glass fiber a bid of 0.2 
mm diameter is used for the present experiments. Care is taken during the 
combustion process such that the droplet remains at the center of the chamber. Fuel 
droplet is generated by using a specially designed fuel delivery system, which 
consists of a hypodermic tube and piston cylinder assembly. Resistance type igniter 
is designed and developed in the lab for igniting the droplet. The igniter is retracted 
just after the ignition by an electro magnet. The igniter and the electro magnet is 
control through an electrically control timer circuit and a relay as shown in Fig. 3.2. 
A 12 V DC power supply is used both for the igniter and electro magnet. The 
droplet diameter change during the burning is visualized by using shadowgraph. A 
backlight is used behind the droplet. The flame shape and the droplet diameter 
change during the combustion are recorded through a microscopic lens of 100 mm 
focal length and a 24 frame per second Hi 8 CCD camera. The droplet bximing is 
recorded as a movie and extracted to the frames using a computer. There it can be 
analyzed frame by frame using image processing technique. The time between 
consecutive frames is constant, which 1/fps. The diameter of the droplet is 
calculated with respect to the support fiber. 


TIMING a)MrONENT5 PECOUfUNG 



Fig. 3.2 : Timer circuit 
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3.3. Experimental Procedures: 

The experimental procedures are discussed in this section: 

i) The fuel is injected to the droplet holder for making a droplet at the tip of 
the droplet holder. 

ii) The shadowgraph system is made operational for the experiments. 

iii) The CCD video camera is focused on the droplet by using the microscopic 
lens. 

iv) The recording of the droplet photograph is initiated. 

v) The ignition system is switched on by triggering the ignition circuit. 

vi) The burning history of the fuel droplet is recorded till its end in video 
format. 

vii) The frames are extracted from the video images by using image processing 
software. 

viii) These photographs of flame and droplet surfaces are modified to the gray 
scale images of 256 steps (0-255) in brightness. 

ix) The flame and droplet boundaries are obtained by examining the contrast 
level of the pixels by using edge detection technique. 

x) The droplet diameter and the flame boundaries are analyzed. 

xi) As at the starting of the combustion there is the igniter, the droplet burning 
is being affected due to this. Hence the initial images have been removed 
and not considered in the diameter regression and flame analysis. 

3.4. Results and Discussions : 

The photograph of flame and the droplet during the combustion of n-heptane with 
initial diameter of 0.82 mm and n-Decane with initial diameter of 1.0 mm are 
shown in Fig. 3.3 and 3.4 respectively. It can be noted that the experiments is 
conducted under normal gravity conditions. Hence, the hot combustion products are 
moved in upward direction due to natural convection resulting in an oval flame 
shape as shown below. The upper part of the flame for both the case of n-Heptane 
and n-Decane burning is very luminous as compared to the lower part. The yellow 
color of the flame can be due to the presence of soot particles. Below this zone there 
is a small zone of reddish color which may arises due to radiation from CO 2 and 
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water vapor. The lower most portion of the flame is observed to be blue due to 
complete combustion in the presence of sufficient oxygen level. As a result the 
maximum temperature is likely to occur in this zone. 



Fig. 3.4: Burning history of n-Decane droplet of diameter 1.0 mm 


Due to buoyancy and the surface tension of the fuel and the support fiber, 
the droplet is not spherical. It is like as shown in Fig.3.5. The droplet diameter is 
approximated by equivalent volume diameter of the droplet as stated by Kumagai 
and et al. [31]. If the major axis diameter is a and the minor axis is b, then the 

equivalent diameter isd = . The schematic representation of the flame envelop 

during droplet burning is shown in Fig. 3.5. The bottom edge of the flame envelop 
from the center of the droplet along its radius is represented by hi as shown in Fig. 
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3.5. Since the droplet flame envelop is non spherical in geometry due to the natural 
convection, another dimension needs to be defined for characterization of the flame 
shape (Kumagai et al. [3 1]). For the present case the side edge of the flame envelop 
from the droplet center along the radius is represented by ha as shown in Fig. 3.5. 
Digitally captured photographs of flame and droplet surfaces were modified to the 
gray scale images of 256 steps (0-255) in brightness. The flame and droplet 
boimdaries are obtained by examining the contrast level of the pixels by using edge 
detection technique. 



h2 

Fig. 3.5: Schematic presentation of droplet burning. 

As per classical theory of droplet combustion the droplet diameter square 
difference, {dl-d^^ is proportional to burning time. The normalized diameter 

square, (f/do^ is plotted against the normalized time, t/do^ in Fig. 3.6 for n-Heptane 
droplet burning of initial droplet diameter 0.82 mm. It can be observed from the 
Fig. 3.6 that droplet diameter decreases with time as per the (£ law of Classical 
theory. The slop of the curve (^/do^ vs. t/do^ known as the burning rate constant, Kb, 
become almost constant in this case which is equal to 1.12 with the experimental 
error is within 4 %. However, the value of Kb obtained from the experimental data 
of Kumagai & Isoda [29] for n-Heptane droplet burning of diameter 1.3 mm at 
normal gravity is approximately equal to 1.04. This small difference may be due to 
use of different initial droplet diameter. For bigger size droplet the radiation loss 
becomes more predominant and is resulted in decrease in the burning rate constant. 
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Fig. 3.6: Variation of normalized droplet diameter square with normalized time for 
burning of n-heptane fuel droplet of diameter 0.82 mm. 



Fig. 3.7: Variation of normalized flame boundaries with normalized time for 
burning of n-heptane fuel droplet of diameter 0.82 mm. 

The flame boundaries, hi/do, Wdo, obtained from digitally captured as 
described in previous section are also plotted in Fig.3.7. It can be observed that the 
flame boundary initially increases with time and subsequently decreases with the 
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droplet diameter. This may be due to fact that tuel gets vaporized in the beginning 
with the initial ignition energy. Subsequently it diffuses outwardly and comes in 
contact with the ambient oxygen at stoichiometric ratio leading to the formation of 
the flame surface at a shorter distance from the surface. As the burning of the 
droplet continues further, flame surface moves outwardly in search of requisite 
oxygen for complete burning. The flame boundaries continue to move outwardly till 
a threshold flame radius beyond which sufficient oxygen is available for complete 
combustion. After this the flame radius starts decreasing with the reduction in 
droplet diameter. 



t/dg (sec/mm^) 

Fig. 3.8: Variation of normalized droplet diameter square with normalized time for 
burning of n-decane fuel droplet of diameter 1 .0 mm. 

Normalized droplet diameter square and the flame boundaries are plotted in 
Fig. 3.8 and 3.9 respectively for combustion of n-Decane fuel droplet of initial 
diameter 1.0 mm. The experimental uncertainty lies within 5 % in this case. The 
burning time for n-Decane fuel droplet is found to be more as compared to n- 
Heptane fuel droplet. This difference in burning time may be attributed to the lower 
volatility of n-Decane. Interestingly the flame of heptane fuel droplet is established 
at a smaller distance from the droplet surface as compared to that of the decane. 
This phenomenon can also be explained in terms of the gasification rate. As it is 
already known the gasification rate is higher in case of heptane fuel due to its low 
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boiling point temperature. The more amount of fuel vapor is likely to be released 
from the droplet surface which spreads for longer distance to get consumed by the 
oxygen. Hence, the flame surface for n-Heptane is located at a larger distance from 
its center as compared to the n-Decane. 
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Fig. 3.9: Variation of normalized flame boundaries with normalized time for 
burning of n-decane fuel droplet of diameter 1.0 mm. 



Fig. 3.10: Variation of normalized droplet diameter square with normalized time for 
burning of 50 % 50% n-heptane n-decane mixture droplet of diameter 0.86 mm. 
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Fig. 3.11: Variation normalized flame boundaries with normalized time for burning 
of 50 % 50% n-heptane n-decane mixture droplet of diameter 0.86 mm. 


An attempt is taken to study the multicomponent fuel droplet combustion. A 
fuel mixture of 50% n-Heptane and 50% n-Decane is taken. The normalized droplet 
diameter square and the flame boundaries for combustion of 50 % -50 % n-heptane 
and n-Decane mixture droplet of diameter 0.86 mm are respectively plotted in Fig. 
3.10-1 1. The results found out here qualitatively similar to the results of Wang et al. 
[76]. The stepped cf law found out in the case of 70% heptane and 30% hexadecane 
may be due to some numerical error. The (f-\dm behavior of burning is observed 
here. The burning rate constant is found to be almost equal to 1.1 with the 
experimental error is within 5%. This is in between the burning rate constant for 
single n-Heptane and n-Decane burning. The droplet boundaries are also in between 
the droplet boundaries found out for these two single pme fuel buniings. 
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3.5. Conclusions 

The experimental studies on single and multicomponent fuel droplet combustion are 
carried out under normal gravity conditions at normal pressure and temperature. 
The droplet diameter regression for both single component and multicomponent 
fuel burning follow the <f law behavior. Flame boundaries are first increase and 
subsequently decreases with time. The slope of the law curve known as burning 
rate constant, Kb, are almost constant for droplet burning. For n-Heptane and n- 
Decane droplet burning under normal gravity the burning rate constant are found to 
be 1.12 and 1.03. In the case of multicomponent fuel, composed with 50% heptane 
and 50% decane, droplet combustion the burning rate constant lies in between the 
burning rate constants of the two fuel droplet combustion. 
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Chapter 4 

MATHEMETICAL MODELING 

4.1. Problem Definition 

In the present chapter, the mathematical modeling of droplet combustion is 
described. The schematic of the model for single isolated fuel droplet is shown in 
Fig. 4.1. In this case, droplet is burning in a hot, stagnant ambience with the 
assumption of spherical symmetry, which is possible only when there is no 
buoyancy force due to gravity. Here the burning is considered to be initiated by 
creating a spherically symmetric high temperature zone around the droplet. Due to 
high temperature around the droplet, the vaporization is started from the surface of 
droplet. Subsequently fuel vapor gets transported by diffusion process and comes in 
contact with the ambient oxygen. Since the gas phase is at high temperature, 
exothermic reactions take place where the fuel vapor and oxidizer meets at their 
stoichiometric ratio. And this zone is defined as flame, which is located at flame 
stand off distance {rj) from the droplet surface as shown in Fig. 4.1. Both the heat 
and the products from this flame zone are transferred to the either sides due to heat 
and mass transfer. The heat transferred from the flame surface causes the droplet to 
vaporize further. For modeling of the droplet combustion one dimension, unsteady 
governing equations such as mass, momentum, energy and species equations are 
considered with the following assumptions. 



Heat & Products 


Flame 
Fuel Droplet 
Fuel Vapor 

Oxidizer 


Fig. 4.1: Schematic of a single isolated liquid droplet burning in an infinite 
oxidizing medium. 
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Assumptions 

i) The droplet is stationary in the quiescent atmosphere. 

ii) The droplet combustion process is considered to be spherically symmetric 
during whole burning period. This assumption is more appropriate only 
when burning takes place under micro-gravity conditions. 

iii) One dimensional unsteady inviscid flow. 

iv) The pressure is assumed to be constant, so that the momentum equations 
need not to be solved. 

v) Dufour-Soret effects and viscous effects are considered to be negligible. 

vi) The combustion is modeled using single step Arrhenius second order 
reaction for simplicity. 

vii) Soot radiation is neglected. 


In the present model variable thermodynamic and transport properties are 
employed to achieve the realistic prediction. A simple gray gas radiation model is 
considered for gas phase radiation. 


4.2. Liquid Phase Conservation Equations 


4.2.1. Energy Equation 


The droplet is heated only by the heat supplied from the flame surface. As the 
interior temperature of the liquid droplet is low, the radiation from the liquid phase 
can be neglected. Liquid phase velocity is also very low as compared to the gas 
phase velocity. Hence it can be neglected. Considering above points the energy 
equation for the liquid phase becomes 


dt 


r^PiCpj 8r\^ dr J 


( 4 . 1 ) 


Where, /?,, T, , are liquid phase density, heat conductivity and specific heat 
respectively. T, is the liquid phase temperature. 
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4.2.2. Liquid Phase Species Equation 


For single component fuel droplet combustion gas phase solubility is negligible at 
low pressure. Hence, there is no need to solve species equation for the liquid phase. 
However, the species equations are to be solved for multicomponent fuel droplet 
combustion. Several models are reported in literature for multicomponent fuel 
droplet combustion. In Diffusion Limit model, the liquid phase mass transfer is 
occurring only due to the molecular diffusion. The governing equation for this 
model is as follows; 


Pi 


^ = 1 
dt 


dr 




( 4 . 2 ) 


where, , is the mass fraction of species in liquid phase. is the liquid phase 

mass diffusivity of f'’ species in liquid mixture. 

During the combustion of multicomponent fuel droplet, the high volatile 
fuel as compared to the other constituent fuels vaporizes at a faster rate from the 
outer shell of the droplet. As a result the scarcity of high volatile substance is likely 
to be occurred at the outer shell. This high volatile fuel vapor can not be transported 
from the inner shell only due to diffusion process. Hence, only the low volatile fuel 
is likely to be present in the outer shell. In contrast both the high and low volatile 
fuels are present in the inner shell. Then there is a chance that high volatile 
component fuel gets superheated due to transfer of heat into the core region of the 
droplet. As a result high pressure can be generated in the core region due to 
vaporization at high pressure, leading to formation of internal circulation inside the 
droplet. This internal fuel circulation can create uniformly distributed species 
concentration over the spatial domain of the liquid droplet. This kind of model is 
known as Distillation Curve Model, which is described as follows; 


^4 3 ^ 

PiJij 


V- 


yy+1 


'St/ 
V 3 


-(Anr^r^M ( 4 . 3 ) 


here, r^ is the droplet radius. £•, is the fractional mass gasification rate for species i , 


which is defined as. 


£i = 


r \ 

% 

V ' = r 


( 4 . 4 ) 
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The first term in the Eqn. 4.3 denotes the total mass of the fuel at (/+1)'^ 
time step, which is equal to the remaining mass of the fuel at the time step after 
vaporization of the fuel from its surface. 


4.3. Gas Phase Conservation Equations 
4.3.1. Mass Conservation 


The mass conservation equation for gas phase is given below; 

■(P.r\) = 0 


dp, , I d ^ 2 


-+- 


dt dr 

where, is the gas phase density, v is the gas phase radial velocity. 


(4.5) 


3.3.2. Species Conservation Equation 


The gas phase species equation is, 
P 


( dY . 

gyt 


dY 

— + v ^ — — 
dt ^ dr j 


r^ dr 


( 


5F„ 




g,‘ 




■ 0 ). 


S<‘ 


(4.6) 


here, Y^ ^ is the species in the gas phase. , is the generation of species due 
to reaction. , is the gas phase diffusional coefficient of species i . 


4.3.3. Energy Conservation Equation 

Considering heat generation due to exothermic reaction, heat conduction, radial 
convection and gas phase radiation of heat, the energy conservation equation for gas 
phase becomes; 


P.C, 




dT^ 

dt 


-+v, 


ar.'i 1 d { 

r^dr\^ dr 


dr 


(4.7) 


Here,C/, g ,'^-gare the specific heat and heat conductivity of the gas phase mixture. 

V-^^are the enthalpy of formation for species and radiation term 
respectively. 

_ 1 5 / 2 X 

r or 


(4.8) 
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4.3.4. Equation of State 


Considering ideal gas, the pressure, temperature, density and mass fraction can be 
related using the equation of state as given below; 



4.4. Reaction Modeling 

During the actual combustion process, several intermediate species are participated 
in the large number of chemical reactions. Some reaction occurs at very first rate, 
some very slowly. The generalized chemical reaction of arbitrary complexity can be 
represented by the following stoichiometric equation given below; 

(4.10) 

where, v] is the stoichiometric coefficient of the reactant, v. is the 

stoichiometric coefficient of the f’’ product, M is the molecular weight of all 
chemical species, and N is the total number of components involved. Then the 
reaction rate, RR can be expressed as 

RR^kfliC^f (4.11) 

/=! 

where, k is the proportionality constant called the specific reaction rate constant. 
For a given chemical reaction, k is independent of the concentrations, and 

depends only on the temperature. The value of k can be determined using Arrhenius 
Law, which is expressed as; 

k = AQxp — ^ (4.12) 

Here, A is pre-exponential factor representing the collision frequency and is the 
activation energy. 


a> = k{Y,X{Yofi 


(4.13) 



33 


In the study of single component fuel droplet combustion, a single step global 
chemistry model for n-Heptane-air is considered as follows; 

C,H,^+\\0^=1C0^+%H,0 (4.14) 

However, for multicomponent fuel droplet combustion, two single step kinetics of 
n-Heptane and n-Decane are used as given below. 


+1 1 02 = 7 CO 2 + 8 H,0 
+ 15.5 O 2 = loco, + 1 1 HP 


(4.15) 


For finding out the reaction rates of the above reactions, the values of A, S, Ea, m 
and n axe taken from Turns [73]. 

4.5. Initial Conditions 

The selections of proper initial conditions are of paramount importance in making 
the numerical process stable. High gradients in the temperature and species profile 
may lead to numerical instability. Hence, suitable initial temperature, oxygen and 
fuel species profiles are chosen (Saitoh and Nagano [61]). Initial profiles with the 
transformed coordinate system are shown in Fig. 4.2 and 4.3. Far away from the 
droplet (at infinity), standard ambient conditions, with standard air composition, has 
been used. At ambient To= 298.15 K with ambient mass fraction values of O 2 , CO 2 , 
HjO taken as 0.232, 0.00045 and 0.0188 respectively. 


At time, t = 0: 


Tgfi = fn(r) , 7] 0 = Mr ) , = /n(r) , = /«('") 

=0.00045, = 0.0188, 


(4.16) 


4.6. Boundary Conditions 
4.6.1. At Droplet Center (r=0) 

Due to symmetry at the droplet centre, the temperature and species concentration 
gradients axe zero. Hence, the boundary conditions can be written as 


(4.17) 
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4.6.2. At Droplet Surface ir=r^ 


Mass Balance: At the droplet surface the liquid fuel gets vaporized into gas 
satisfying the mass conservation as stated below; 


=-Pr 


(4.18) 


Species Balance: In both of the phases, species are transported due to both 
convection and diffusion. Therefore, species flux is governed by; 


‘ *■' dr 


— p V Y 


/ = 02,C02,H20 (4.19) 


= Pg'^g(^f-^gJ,r-r,) f = ^f2 


(4.20) 


-PAj- 




(4.21) 


Energy Balance: At droplet surface certain amount of heat from the gas phase is 
transported to the liquid droplet. The heat balance equation at the liquid-gas phase 


interface can be written as. 


dT dT, dr. A 

K — ~ = ^1 — ~ - Pi ~ / A + 

^ dr ‘ dr dt ^ ‘ ‘ 


(4.22) 


The first term is pertaining to heat transfer due gas phase conduction; 
second term represents heat absorption by the liquid due to heat conduction. While 
the third term represents the latent heat of vaporization. The forth term is associated 
with the radiation. 


Phase Equilibrium at Liquid Vapor Interface: The phase equilibrium at the 
liquid vapor interface of the droplet surface is governed by the Clausius 
Clapeyron’s equation and the Roult’s Law as stated below [54]; 


1 LMA 1 1 

P A/ I ■‘A./ h 


(4.23) 


Here, A is the surface temperature, A ^ is the boiling point temperature 

evaluated at the normal atmospheric pressure, P is the prevailing system pressure in 
bar and and fyare the molecular weight and latent heat of vaporization 
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respectively for the i species, and is the universal gas constant. It can be noted 

that the above expression for phase equilibrium at liquid vapor interface is only 

valid for ideal gas case. Hence, it cannot be used for combustion process at high 
pressure. 

4.6.3. At Infinity (r=oo) 

At the ambient, the temperature or species profiles are not affected by the droplet 
burning. Hence at the ambient, the temperature and species profile remain constant 
throughout the burning process, which is expressed as; 

Y^j=0 (4.24) 


4.7. Non-dimensionalization and Co-ordinate Transformation 


4.7.1. Non-dimensionalization 


The droplet surface regresses with the burning of fuel droplet, which creates 
additional complexities in the modeling. In order to overcome this complexities co- 
ordinate trans-formation is required. The governing equations are non-dimension- 
alized by using the following non-dimensional terms; 







( Cp ufiCj- 
\ ^co J 




• r\ ^ 

V 0= = — 


'ip N 

i=\ 


Y:=a,Y, (4.25) 


where, the co-efficient, a, is known as absorptivity and can be defined as follows: 

M N N 

a and M = (4.26) 

(V,. -V,.)M,. 

All thermo-physical properties are non-dimensionalized by their reference 
value taken at the ambient condition. The non-dimensionalized expressions are; 
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The ajand D^ Qare further non-dimensionalized with the Lewis No. of the 
species i, le, 


(4.28) 

The asterisk (*) mark indicates the non-dimensional quantities. 


4.7.2. Co-ordinate Transformation 


In order to overcome the problem of moving boundary at the droplet surface it is 
convenient to transform the moving boundary to a fixed boundary. In liquid phase 
the gradients do not vary much over the spatial domain, however, for gas phase the 
gradients are found to be quite high between the droplet surface and the flame. 
Hence, to compute efficiently, very fine grids are to be provided within this zone. 
On the other hand the gradients of dependent variables are found to be low between 
the flame and the ambient. Hence, course grid should be used to save computational 
time. To overcome these difficulties the co-ordinate transformations are adopted 


from the work of Saitoh and Nagano [61]. 

For liquid phase the following transformation is employed; 


(4.29) 


Note that, the above transformation converts the liquid region to a fixed 
domain from x = 0 (at the droplet centre) to x = 1 (at the droplet surface). 

For the gas phase, the following transformation is used. 



This transformation converts the gas phase to a fixed domain. At the droplet 
surface, 17 = 0 and 7 = 5 is considered at far distance {r ->oo) from the center of the 
droplet. With this transformation the grid become finer between droplet surface and 
the flame and coarser beyond the flame region. 
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4.7.3. Transformed and Non-dimensionalized Governing Equations 

i) Non dimensional energy equation in liquid phase: 

The energy conservation equation (Eq. 4.1) is transformed to 


ar; _ a' dX 
dt* r? cbc' 


+ 


1 


ai; 


2a] 

xr/ ' p] C] j r/ dx 


dT' 


dx 


(4.31) 


ii) Non dimensional species equation in liquid phase : 

The diffusion limited species conservation equation (Eq. 4.2) is transformed 
to 


dYj Dl d% ^ r 2Dl ^ 1 dP], 
dt* LeXd^ ^ ^ 


X. 


dx 


(4.32) 


The equation of species profiles for distillation curve (Eq. 4.3) model 
becomes 

/ Ptd K> 1 pIj 1 - ( [^d )' ^d pI ) (4-33) 

Jj+\ V-3 Jj ' 'J 


Hi) Non dimensional Gas phase mass conservation equation: 
The mass conservation equation (Eq. 4.5) becomes, 


dt 




'd j 


dq r 


dq 


(4.34) 


iv) Non dimensional gas phase species conservation equation: 
The species conservation equation (Eq. 4.6) can be derived as 

1 




dt 


Le^or^ drj 


+ 




+ 


Le,X" r* Le,X"pl dq 




dq 


(r;)"(r4)"exp 


T* 

V ''s 7 


(4.35) 


v) Non dimensional gas phase energy conservation equation: 

The transformed non dimensional energy conservation equation (Eq. 4.7) is 
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dj,_ < dX 


dt 
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dt]^ 


+ 
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3W )' 


dr, 


dT' 
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*3 

^ P. 


^ ^ / *-2 V Dn / frr* 

;r^— (^ ?«)+-^(p,) {Yf) (yX exp 

g ^P,g ^P,g 


m+n-\ 


0 


V 


(4.36) 


where, the Damkohler number, Dn , is defined by, 

[^fMfJ(cCoMo,)" 


(4.37) 


4.7.4. Transformed and Non-dimensionaiized Initial Conditions 

For liquid phase at t* = 0, 

t;=t;{x), y;,^y;j^, for o<x<i 

For gas phase at t* = 0, 

Kj=KA^)p Ko.=K,oSv) for0<7<co 


(4.38) 

(4.39) 


Initial temperature, fuel and oxygen species are chosen as shown in Fig.4.2-3 with 
the transformed coordinate system grids. 



Fig.4.2: Initial Temperature Profile for (a) Liquid Phase and (b) Gas Phase. 
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Fig.4.3: Initial Oxygen and Fuel Concentration Profiles. 


4.7.5. Transformed and Non-dimensionalized Boundary Conditions 

(i) At droplet center, jc = 0 , the equation 4.17 becomes, 

at, 


dx 

dT; _ 

dx 


= 0 
= 0 


(4.40) 


(ii) At droplet Surface, jc = 1 , mass conservation equation (Eq. 4. 1 8) is converted 
to 


♦ • * • * dft 


(4.41) 


(iii) At droplet surface species conservation equations (Eqs. 4.19-21) can be 
transformed to; 


dYl 




dx 


L^i 0 

— — v„7„,. 


g -g,/,r=o 

x = l «■' 


i = 02,C02,H20 (4.42) 


dY 


zJ 


dx 


K. 


Le^n r 


x = \ 


D 






dx 


x = l 




f = fvf2 


(4.43) 


(4.44) 
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(iv) Energy balance equation at the droplet surface (Eq. 4.22) becomes 


5 ^ ^ St 

dt’ r 


‘ dx ^ dr] 




+ <1 r 


(4.45) 


where, St is a non-dimensional number defined as: 

St — ^co 




(4.46) 


/=! 


4,8. Numerical Technique (Discretization) 

One dimensional, unsteady, equations for mass, species and energy equations for 
both liquid and gas phases governing the combustion processes of a single droplet 
fuel are solved numerically as described below. The Cranck-Nicolson finite 
difference scheme (Anderson, [2]) is employed for the discretization of governing 
equations. The unsteady term can be discretized as follows; 


dT 



rj^n 


dt At 


(4.47) 


The second derivative terms are discretized by using central time step, 
central spacing (CTCS) as described below; 

d^T _ 1 
dt]^ 2 

Diffusion terms are also expressed by using central time step, central spacing. The 
discretized diffusion terms is shown below; 

^ = 1 
drj 2 


r ^« + l 

■‘(+1 -‘>-1 


V 


lArj 


. + 


‘/+i 


-T 

■‘(-1 


\n 


) 


2Ar] 


(4.49) 


/ \«+l 

r _ 

T,.,-2T, + T„ 

+ 



1 J 

1 J 



(4.48) 


The terms, — , — are discretized by using backward time spacing, for the 
’ dt] drj 


^3 -3 

A+1 A-i 


2 At; 


previous time steps, as; 


dt] 


V 


(4.50) 
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For discretizing of the species equation at the droplet boundary, the 
following discretizations are employed. 


dr] 


77 = 0 


/^y_y^« + l ^ 
■‘2 -'1 


A 7 


+ 


V 


y 


(4.51) 


For solving the mass conservation equation for the gas phase the interior 
grid points are discretized by higher order central spacing and the points adjacent to 
boundary are discretized using second order central spacing technique. 

dp _ -3/?, +4/2,„ - A., 
dr] 2 A 7 


_ ~Pm + P/-2 

drj \1^T} 

67 2 A 7 


i = 2 

/or i = 3to m-2 (4.52) 

for i-m-l 


The non dimensional energy conservation equation for gas phase is 
discretized using Crank-Nicolson finite difference scheme as given below. 
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J for 2 = 2torrt-l (4.53) 


-r + <p 


In the above equation is replaced by T for convenience. The terms 
a,p,y and cp are defined as follows; 
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Equation 4.53 becomes 
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(4.58) 


2A77‘ 4At7 


7’”, -^Ar* +^A/’ 


In these way, this equation can be represented as a tri diagonal matrix form 
as shown below. 
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a At 
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2Ar] 4A7 


Tl,-rAt*-r(pAt* (4.60) 


The values of a,, 6,, a„,c„ can be obtained from the boundary conditions. 

This matrix is solved very efficiently by using Thoma 5 Algorithm. Similarly the 
other governing equations for species and energy, (Eqs. 4.31, 4.32, 4.35) with 
appropriate boundary conditions are resulted in tri-diagonal matrices and are solved 
by using the same algorithm. The solution procedure for the present model is 
described in details in section 4.1 1. 


4.9. Radiation Heat Transfer Modeling 

For gas radiation, the gray gas model has been used. In this model, carbon dioxide 
and water vapor are assumed to be the major contributors for the radiative heat 
transfer. The Planck mean absorption coefficient kp^is evaluated locally from 
mass- fraction weighted summation of each participating species as folio w's: 


(4.61) 
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Here, the mole fractions of CO 2 and H 2 O, respectively. P is the 

atmospheric pressure in bar. The Planck mean absorption coefficient for CO 2 and 
H 2 O are represented by and respectively. These are determined by the 
curve-fitted data taken from RADCAL program 1993[19]. 

The radiation term involved in the energy equation has been calculated by 
using radiative transfer equation with the assumption of isotropic scattering, given 
by Viskanta and Merriam [75] for a concentric sphere. The same modeling has also 
been used by Saitoh et al., [62], Marchese et al. [43] and Kumar et al., [32], [33]. 




(4.62) 


where, is radiation flux, is the Stephan-Boltzman constant, G is the radiation 

absorption by the gas shell. The first term is associated with radiation emission and 
the second term corresponds to the absorption; however, the scattering term gets 
canceled out. 

The radiation absorption can be expressed as. 


G(t) = ^ 




(4.63) 


Here, g, (r) is associated with the radiation absorption by any shell from the 
droplet surface. The term, corresponds to the absorption of radiation heat 
from the ambience. The Kernel term, K(r,^) represents the radiation absorption 
between the shells, r is the optical path length defined as: 

(4.64) 

0 

r, and axe the optical path length at the drop surface and infinity respectively. 
The functions g,(T) and g 2 (f) axe defined as: 

= (4.65) 

g,(x) = X, E^Cx, - x)-E,(x 2 -x)-x; E,(x; +x)-E,(x; +r*) 


(4.66) 
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Here r* and rj axe as follows; 

, r; = ^(T/-rf) (4.67) 

Kernel K(r,^) is represented by: 

ii:(T,f) = £,(|T-f|)-£,(|r-r|)xf (4.68) 

where, ^ -r,^ J and the exponential integral function, £.(r) is defined as; 

00 

£.(f)=Jr'exp(-rO<if (4.69) 

0 

The radiation properties have been calculated with high accuracy in this 
way. The close form solutions [85] for Exponential Integral Fmictions are used. 

4.10. Thermodynamic, Transport and Radiation Properties 

The gas phase and liquid phase transport, thermal and radiation properties are 
discussed in details in Appendix A. The gas mixture consists of fuel vapor, O 2 , 
CO 2 , H 2 O and N 2 . The transport properties of the liquid and gas phase have been 
estimated from the kinetic theory of gases and liquids as discussed in Reid et al. 
[59] and Hirschfelder et al. [22]. These estimated are compared with the 
experimental data (Vargaftik, [74] and Toulokian, [72]). The species transport and 
thermal data are taken mostly from Sandia Database. 

4.11. Solution Procedure 

The flow chart of the solution is shown in Fig.4.4. For the present study a computer 
code has been developed by using C with this algorithm. 

( 1 ) At r* = 0 , suitable initial temperature and species profiles are chosen. 

(2) The vaporization rate data for initial conditions is estimated from the 
classical theory of droplet vaporization and combustion. 

(3 ) Calculation of variable gas and liquid phases’ properties with temperature. 

(4) Temperature of a single component fuel droplet at its surface can be 
calculated by using the liquid vapour phase equilibrium equations (Clausius- 
Clapeyron’s equation and Roult’s Law). In the case of Multi-component fuel 
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(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 


( 11 ) 


( 12 ) 


(13) 

(14) 

(15) 


droplet surface temperature can be calculated by mass fraction weighted 
summation. 

To find out the axi-symmetric velocity, , mass conservation equation has 
to be solved. 

Determination of the Radiation heat transfer and irradiance. 

Solution of Gas phase energy conservation equation to find out the 
temperature profile. 

Then species equation (fuel, O 2 , CO 2 , H 2 O) has to be calculated respectively 
to find out the species profiles. The species profile of N 2 is determined from 
species balance. 

Calculation of the energy equation for the liquid phase to find out the 
temperature inside the droplet. 

For multi-component fuel droplet combustion the species equation in the 
liquid phase are solved. For distillation curve (rapid mixing) model uniform 
species concentration are considered inside the droplet and for the diffusion 
limit model diffusion limited governing equation for the liquid phase species 
equation have to be solved. 


Determination of the regression rate 




at surface of the fuel droplet 


using the energy conservation equation. 


Check convergence criterion. If 


^ dr ^ 




dt 


ydtj 


<10^ , then go to 


old 


the next step, else, do the same process 3-11. 


Droplet radius for the next time step has to be calculated, rj+, = /j + 


\dtj 


At. 


Check the stopping criterion. 
Proceed to the next time step. 


All the computations are carried out in the present thesis using Linux based 
Pentium IV (3.4 GHz) IBM PC. The time required for a typical case in this 
machine is around 1 CPU hours. 
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Fig. 4.4: Flow Chart of the Solution Procedure 
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4.12. Selection of Grid Sizes, Time Step and Domain Boundary 

Ideally, the outermost boundary lies at infinity. However for numerical simulation, 
outer boundary is generally taken at considerable distance from the droplet centre. 
In the present numerical model the outermost boundary is taken at approximately 
150 times away from the droplet radius (7 = 5.0). In order to choose the 
appropriate semi-infinite domain for computation, the effect of outer boundary 
distances on the results of regression rate are investigated by varying the value of 
7 “4, 5, 6 , for an n-heptane droplet at normal ambient temperature and pressure. 
The results of d^-law for these cases are plotted in Fig. 4.5. It can be observed that 
all three curves are pretty closer to each other. The errors are within 0.01% as 
compared to the total burning time. Hence, the outer boundary disfrince, 77 of 5 is 
chosen for all the computations reported in this thesis. 



Fig. 4.5: The results of law for n-Heptane burning in normal atmosphere by 
varying the boundary distance, 7 = 4, 5 , 6 ( At* = 5 x 1 0"^ and A;; = 0.02 ). 

A uniform non dimensional time step (At*) of 5x10-4 has been used. To 
check the accuracy of time step chosen the following test has been performed. In 
Figure 4.6(a), the / with time for n-Heptane burning in normal atmosphere is 
plotted with, At*, of 2x10“^, SxlO"^ and 1x10"^ with A;/ = 0.02 and at ambient 
7 is considered equal to 5.0. From this figure, it has been fond that At* =5x10'^ is 
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the optimum grid size, as there are negligible differences in the result for 
At -5x10 and At =2xl0~^ (0.2% of error as compared to the total burning 
time). Hence, At =5x10'' is taken for the entire calculation. 



Fig. 4.6: The results of cf law for n-Heptane burning in normal atmosphere with 
(a) Different At of 2x10”'', 5x10"^ and 1x10”^ (Atj = 0.02, at ambient t] is 
considered equal to 5.0) and (b) different Arjof 0.01, 0.02 and 0.04 {At =5x10"'' , 
at ambient rj is considered equal to 5 . 0 ). 

The logarithmic grid is used A 7 = 0.02 . In order to select optimum grid size 
grid sensitivity test has performed. The results of the variation of burning history 
with the different grid size taken are plotted in Fig. 4.6 (b). It is depicted from the 
figure that, burning time does not depend much on the grid chosen. For At]= 0.02 
and 0.04 the relative in within 0.1 %. Hence to save the computational times the 
value of A 77 is chosen 0.02 for entire calculations in the thesis. 
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Chapter 5 

RESULTS AND DISCUSSION 

In this chapter, the results achieved from the numerical studies of single component 
and multicomponent fuel droplet combustion are discussed. The results are 
compared with the present experimental data and others data available in the 
literature. First section (Sec. 5.1) contains the studies of single component fuel 
droplet combustion followed by multicomponent fuel droplet combustion in the 
next section (Sec. 5.2). 

5.1. Single Component Fuel Droplet Combustion 

This section is mainly devoted to the numerical studies of single component fuel 
droplet combustion. In the first phase, the numerical experiments are carried out 
without including the radiation effects. A single n-Heptane fuel droplet of initial 
diameter of 1 mm is being investigated. The results of this study is presented in 
subsequent subsection in terms of droplet burning history, flame stmcture, flame 
stand-off ratio and flame location and gas and liquid phase temperature distribution. 
The numerical data such as droplet burning history and the flame location are 
compared with the present experimental data and other experimental data available 
in the literature. In the next subsection the effects of radiation heat transfer are 
being discussed extensively. The effects of initial droplet diameter on the burning 
rate at different ambient temperature are investigated and compared with the 
existing experimental data. Subsequently the results of droplet burning at different 
ambient oxygen level are presented. 

5.1.1. Transient Variation of Burning Parameters without Radiation 
Droplet Burning History 

The normalized droplet diameter square, ^/do and the burning rate constant, Kb 
are plotted with normalized time, t/do in Fig. 5.1(a) and 5.1(b) respectively. The 
variation of droplet diameter square with time is compared with the numerical data 
presented by Saito & Nagano [61] and the experimental data by Okajima and 
Kumagai [52]. It can be observed from the Fig. 5.1(a) that, the present numerical 
prediction burning rate is higher than the experimentally observed bummg rate by 
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Okajima and Kumagai [52]. In contrast, the numerical results of burning rate 
considering constant property and without radiation predicted by Saito and Nagano 
[61] matches well with the experimental data of Okajima and Kumagai [52]. In the 
present experimental study under normal gravity the burning rate is observed to be 
higher than the microgravity experimental data. However, the present numerical 
study without radiation heat transfer is closer to the present experimental data. 

It can be noted that, radiation heat loss is not considered in the present 
model. Hence, burning rate must be greater than the experimental results, as 
experimental data can not be free from radiation heat loss. In contrast the numerical 
results presented by Saito and Nagano [61] predict well the burning rate 
experimental data of Okajima and Kumagai [52]. This is contrary to the 
expectation. This matching of numerical prediction with the experimental data is 
considered to matter of mere coincidence. However, the present model using 
variable properties are supposed to be superior as compared to the numerical 
modeling with constant properties of Saito and Nagano [61]. 

In the present experiments under normal gravity, natural convection is 
occurred due to which turbulence is likely to be generated in the gas phase. Hence, 
gas phase heat and mass transfer are increased to a larger extent resulting in 
increase in droplet burning rate. The increase in burning rate under normal gravity 
condition is mainly contributed by the natural convection rather than the radiative 
heat loss. 



Fig. 5.1: Variation of the (a) Droplet diameter square and (b) Burning rate constant, 
Kb with time. 
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The transient variations of droplet center temperature, droplet surface 
temperature and flame temperature are plotted in Fig. 5.2. It can be noted that the 
numerical calculation of droplet combustion are initiated with some initial 
temperature and species profile. 



Fig. 5.2: Variation of droplet surface temperature, droplet center temperature and 
flame temperature with time. 


Fig. 5.2 indicates that the temperature profiles change its slope indicating 
the effect of initial transients. The heat generated through the combustion is used to 
build up the temperature in the gas phase. However in this period the flame 
temperature as well as the droplet surface temperature are increased to a certain 
value and become almost stable. During this period the burning rate constant, Kb, is 
increased at a faster rate and attains a maximum value as can be seen in Fig. 5.1(b). 
This period is defined as initial preheated period. In the numerical study this initial 
preheat period is a depended on the initial temperature and species profiles. After 
initial transient period, the droplet burning rate, Kb remains almost constant. After 
this, the gas phase temperature increases marginally to its asymptotic value as can 
be seen in Fig. 5.2. During this period the heat generated due to reaction is almost 
balanced by the heat transferred to the ambient and liquid phase. As a result the 
burning rate constant, Kb, remains almost constant. The droplet surface temperature 
as shown in Fig. 5.2 becomes almost constant during this period. However the 
dronlet center temnerature i?. increased almost at a constant rate as nearly constant 
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amount of heat is being supplied to the liquid phase from the gas phase side. This 
period of burning can be defined as equilibrium burning period or J^-law period. At 
the end of the burning, heat generated by the gas phase combustion may not be 
sufficient to meet the demand of heat losses for the self sustaining of the 
combustion wave, which leads to the extinction of the flame. This period is defined 
as flame extinction period. As a result the burning rate constant, Kb, is decreased 
drastically as shown in Fig. 5.1 (b). Based on the above discussion, the droplet 
burning regime can be divided into three major periods such as (I) Initial pre-heated 
period, (11) Stable <i^-law or equilibrium period and (III) Flame extinction period. 

Flame Structure 

Typical structure of flame at {t/d^) = 0.4 sec/mm^ is shown in Fig. 5.3(a). It can be 
observed in this figure that the fuel and oxidizer species meet at the point where 
maximum temperature occurs at which flame is located. In this zone exothermic 
reactions take place as shown in Fig. 5.3 (b). As shown in the Fig. 5.3(a)-(b), very 
near to the droplet surface {ln{r/r^ towards zero), fuel concentration is higher and 
oxygen concentration is zero. In this zone fuel vaporizes from the droplet and 
subsequently diffuses to the gas phase. However, oxygen is also diffiised to the gas 
phase from the ambient. At a certain distance from the droplet surface fuel and 
oxygen meet at their stoichiometric ratio. In this zone heat is being generated 
resulting high temperature. This high temperature zone, where fuel vapor and 
oxygen meet at their stoichiometric ratio, is defined as flame. The heat is being 
transferred to both the sides from the flame. Due to the reaction between fuel and 
oxidizer product species are generated at the flame zone. These product species are 
convected and diffuses to the ambient as well as towards the droplet surface. The 
typical variations of concentration of the product species are shown in Fig. 5.3(a). 
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Fig. 5.3: (a) Structure of flame and (b) Radial distribution of reaction rate at the 
middle of droplet burning {t/do^ = 0.4 sec/mm^) 

Flame Standoff Ratio and Flame Location 

The variation of predicted flame stand off ratio during the droplet burning is plotted 
in Fig. 5.4(a). The flame standoff ratio increases monotonically with the time and 
then attains a maximum value in the stable law period as shown in Fig. 5.4(a). 
The present predicted flame stand-off ratio data are found to be in the range of 1 1 - 
28. In contrast the numerical data predicted by Saito and Nagano [61] shows the 
flame stand off ratio increases from 4-35. It increases at a lower rate in the initial 
portion of burning but at a faster rate towards the end of the burning. These features 
may be due to the use of constant properties assumptions in their modeling. 

The flame diameter first increases with time and after some time it decreases 
with the droplet diameter. The flame diameter observed from the present simulation 
matches qualitatively with the data presented by Okajima and Kumagai [52]. 
However, the present result fails to predict the experimental data of flame diameter 
quantitatively. This discrepancy might be caused due to the radiation heat loss. 
Although the flame diameter data predicted by Saito and Nagano [61] matches well 
in the initial phases with the experimental data but it departs drastically as can be 
seen in Fig. 5.4(b). It seems that the flame diameter prediction towards the ends of 
burning is quite unexpected as it is known that flame diameter drops down towards 
the end of the burning. This may be due to constant properties assumptions in their 
formulation. 


Nondimensional Reaction rate 
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I^ig. 5.4: Variation of (a) Flame Stand-off ratio and (b) Flame diameter during 
droplet burning. 


CSas Phase and Liquid Phase Temperature Distribution 

fransient variation of gas phase temperature distribution along the radial direction 
is plotted in Fig. 5.5 (a). It can be observed from Fig. 5.5(a) that the gas phase 
temperature increases from the surface temperature of the liquid droplet and attains 
a maximum value at the flame zone. Initially the gas phase temperature increases 
with time at a faster rate. After the initial pre-heated transient period it does not 
vary further with time. The reason may be attributed to the initial heat generation 
due to exothermic reactions, which is used to build up the temperature in the gas 
phase. After initial transient period it is expected that the heat generation due to 
chemical reactions is balanced by the heat loss to the droplet surface and ambient. 
In F'ig. 5.5(b), the transient variation of liquid phase temperature is plotted along the 
radial direction. At the interface of liquid droplet and gas phase the surface 
temperature increases with time initially and attains a constant value. It can be 
observed Irom the Fig. 5.5(b) that the temperature gradient along the radial 
direction with time and attains almost uniform temperature profile across the liquid 
droplet towards the end of burning. This feature is well predicted as per the heat 
transfer due to conduction. 
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Radial Distance, n = (M>-/r;) Radial distance, a: = (r/r; 


Fig. 5.5: (a) Gas phase and (b) Liquid Phase temperature distribution 


5.1.2. Effect of Ambient Oxygen Concentration on Droplet Burning 

The oxygen level in the actual combustion zone is likely to be changed with 
time. Hence it is very essential to study the effects of the ambient oxygen 
concentration on the droplet combustion. An attempt is made to investigate 
numerically the effect of ambient oxygen level on single isolated n-heptane fuel 
droplet burning. 



t/dj (sec/mm^) (secW) 

Fig: 5.6: Variation of (a) cf law and (b) burning rate constant, Kb with different O 2 
concentration at ambient air. 
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The effect of ambient oxygen concentration level on the droplet burning 
history is plotted in Fig. 5.6(a). It can be observed from the Fig. 5.6 (a) that the 
burning rate increases with the ambient oxygen level. This might be attributed to 
increase in the heat generation rate, which is the function of in the 

single step chemistry model. Similar phenomenon can also be predicted by the 
classical quasi-steady theory. 

It can be observed from the Fig. 5.6(a) that the initial preheat period is also 
affected due to the oxygen concentration level. In this regard, total burning time and 
initial preheat time are plotted in Fig. 5.7(a) and 5.7 (b) respectively with the 
relative oxygen concentration level, which is defined as the ratio of the oxygen 
concentration level at the atmosphere to the oxygen concentration level at standard 
condition. The burning time and preheat time decays exponentially with oxygen 
concentration level as shown Fig. 5.7(a) and (b). Besides this, the total burning time 
is more affected by the oxygen level than the initial preheat time. 



Fig: 5.7: Effect of relative oxygen concentration (O 2 Cone. / O 2 Cone, at standard 
atm) at the ambient on the (a) burning time (t*) (b) Preheat time (t,) 

The transient variation of flame stand off ratio is plotted in Fig. 5.8(a). From 
Fig. 5. 8 (a) it can be observed that, the flame stand off ratio is increased with the 
ambient oxygen concentration level. It is known that the flame stand off ratio is 
controlled by balance between the gas phase radial velocity and the diffusional 
velocity of oxygen. In order to understand the increase of flame standoff ratio with 
oxygen level, the transient variation of the gas phase velocity, v^ , at the droplet 
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surface is plotted in Fig.5.9 (b) for different ambient oxygen concentration level. It 
can be observed that the flame temperature is increased with the ambient oxygen 
concentration level as shown in Fig. 5.9 (a) leading to the increase in the droplet 
regression rate as well as the gas phase radial velocity. Hence, flame is moved 
away further from the droplet surface with increase in gas phase velocity. In 
contrast, the diffusional velocity of oxygen, acting in inward radial direction is 
increased with the ambient oxygen concentration level. It attempts to push the flame 
front towards the droplet surface. However the magnitude of the diffusional 
velocity is less compared to the gas phase radial velocity at the flame surface. 
Hence, the flame stand off ratio increases with the ambient oxygen concentration 
level. Interesting to note that present prediction of flame standoff ratio is in contrast 
with the classical theory. These anomalies may be due to fact that the classical 
theory does not take account of the gas phase velocity, which is the determining 
factor in predicting the flame standoff ratio. 



t/dg (sec/mm^) t/d^ (sec/mm^) 


Fig: 5.8: Variation of (a) flame Stand-off ratio and (b) flame stand-off distance 
with different oxygen concentration of atmosphere. 
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Fig. 5.10. It can be observed that present numerical results matches w^ell within the 
error band of experimental data of Mikami et al. [46]. 



Fig.5.10: Variation of flame temperature as a fimction of time {t/do)'- predicting 
along with experimental values (Mikami et al., [46]) 

The normlalized droplet diameter square is plotted with normalized time 
{t/do ) in Fig. 5.11(a). It is found from the Fig. 5.11(a) that the burning time is 
increased considerably due to the inclusion of radiation model. The data presented 
in the Fig. 5.1 1(a) are also compared with the compiled experimental data presented 
by Kumagai et al. [31] and the numerical prediction by Marchese et al. [46]. Only 
the stable law portions data are compared. Interesting to note that, by the 
inclusion of radiation model the present results of burning rate predicts well with 
the experimental data of Kumagai et al. [31] and numerical data of Marchese et al. 
[46] as shown in Fig. 5.1 1 (a). 

The variation of burning rate constemt, Kb at ambient temperature is plotted 
in Fig. 5.1 1(b). It can be observed that by inclusion of radiation model the burning 
rate constant, Kb decreases substantially as compared to results without radiation 
model. However when the ambient temperature is increased to 1093 K, the burning 
rate constant increases substantially as shown in Fig. 5.11 (b). Similar feature can 
also be observed in Fig. 5.11 (a) indicating that increase in ambient temperature 
reduces the burning time. This phenomenon can be attributed to the increase in 
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adiabatic flame temperature with increase in ambient temperature. Beside this the 
heat loss from the flame surface will be hindered due to the reduced temperature 
gradient between the flame surface and its surroundings. 
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Fig. 5.11: Effect of radiation heat transfer on the (a) cf^-Law and (b) burning rate 
constant, Kb at different ambient temperatures for isolated droplet of initial 
Diameter 1.0 mm. 



Fig. 5.12: Effect of radiation heat transfer on the (a) flame temperature and (b) 
flame stand off distances at different ambient temperatures for isolated droplet of 
initial diameter 1 .0 mm. 

An important feature can also be observed from the Fig. 5.11 (b) that the 
burning rate constant increases with time during stable portion of the droplet 
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combustion. However, it remains almost constant for no radiation model. This 
phenomenon can be explained in terms of the effects of geometry on droplet 
burning. The radiation heat loss from the flame front to its surrounding increases for 
larger size droplet due to increase in its surface area. In contrast radiation loss from 
the smaller droplet becomes less. Hence, the radiation heat loss from the flame 
surface is likely to be decreased with reduction in droplet diameter. As a result, 
more amount of heat can be transferred to the droplet surface leading to an increase 
in burning rate. 

The flame , temperature for n-Heptane fuel droplet combustion of initial 
diameter of 1 mm at two different ambient temperatures considering radiation heat 
transferred is compared in Fig.5.12 (a) with the predicted data of present numerical 
model without radiation. A considerable amount of flame temperature drop from 
2200 K to 1700 K can be observed with the inclusion of the radiation model. By 
increasing the ambient temperature to 1093 K the present model predicts the flame 
temperature to be around 2800 K. But, this predicted flame temperature seems to be 
quite high for which gas phase dissociation reaction model may need to be 
considered. 


2.5x10* 
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Fig. 5.13: Transient variation of radial distribution of (a) radiation heat loss and (b) 
heat gain in the gas phase. 

In order to access the mechanism of radiation heat exchange during the 
droplet combustion, the variations of radiation heat loss and heat gain along the 
radial direction are plotted in the Fig. 5.13 (a) and 5.13 (b) respectively at different 
time intervals. It can be observed from the Fig. 5.13 (a)-(b) that, the maximum 
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radiation heat loss is occurred at the flame zone due to high temperature and it 
varies with the temperature. However the more radiation heat gain is likely to be 
occurring towards the droplet surface. Radiation heat gain is maximum at the 
droplet surface. It also increases with the time as the droplet diameter decreases. 


Influence of Initial Droplet Diameter at Different Ambient Temperature 

In order to study the influence of initial diameter on the n-Heptane fuel droplet 
combustion, three different representative initial droplet diameters such as 1.0 mm, 
1.5 mm and 2.0 mm are considered at two different ambient temperatures. Figures 
5.14 and 5.15 show the typical variation of normalized droplet diameter square and 
the burning rate constant, Kb, with time at two different ambient temperature of 
298.15 and 1093 K respectively. The results are compared with the experimental 
data of Xu et al. [79], 


<N o 


1.0 

» 1 '"1 ' ' ’ I 

Initial Diameter ' ■ 
1 .0 mm 

0.8 


1.5 mm 

2.0 mm 

0.6 

\ 

\ ’’n 

\ m 

■ 0.94mm[79] 

V L46mm[79]- 

0.4 

\ 

: n-Heptane 


0.2 

■T = 298.15 

amb 

m 

\ '57 

\ . 

0 n 

^ (a) 

,1 1, 1 ,1 ,1 1 

I ... 1 . . . . ■ 1 i i ., ■ 1-.1: 


t/d^ (sec/mm^) 



Fig.5.14: (a) cf law and (b) Burning rate constant for different droplet size (1.0 
mm, 1 .5 mm and 2.0 mm) at room temperature 


It can be observed from the Fig. 5.14 that the burning rate decreases with the 
droplet diameter at low ambient temperature, 298,15 K. As a result the burning 
duration for complete combustion of droplet increases with the increase in droplet 
diameter. As the droplet diameter is increased the flame is located at a larger 
distance from the surface of the droplet. As a result, for a bigger droplet the heat 
loss from the flame surface to the ambient is likely to be more as compared to the 
smaller size droplet. Hence this is resulted in decrease in burning rate constant. 
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Fig.5.15: (a) ^ law and (b) Burning rate constant for different droplet size (1.0 
mm, 1 .5 mm and 2.0 mm) at ambient temperature = 1093 K. 

In contrast the predicted variation of burning rate at higher ambient 
temperature of 1 093 K exhibits an opposite trend with increase in droplet diameter 
in comparison to the low ambient temperature case as shown in Fig. 5.15. Because 
at high temperature the heat gain from its surrounding is expected to be more as 
compared to the heat loss. This effect is more pronounced for bigger size droplet. 
Hence, the initial diameter effect has an opposite trend at high ambient temperature. 
Similar features can also be observed in the experimental data of Xu et al. [79]. 
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Fig. 5.16: Influences of initial droplet diameter on the average burning rate constant 
at different ambient temperature. 
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In order to compare the predicted burning rate constant with the 
experimental data available in the literature, the burning rate constant is averaged 
over the stable burning period and is plotted in Fig. 5.16 against the initial droplet 
diameter at two different ambient temperatures (298.15 K and 1093 K). The 
experimental data of Okajima Sc Kumagai [52], Kumagai et al. [31] and Xu et al. 
[79] are plotted in this figure. It can be observed that the present prediction matches 
well with the experimental data of Kumagai within the experimental error band. 
Although the experimental data of Xu et al. [79] exhibit the similar trend than that 
of the present numerical data and experimental data of Kumagai, but it is in the 
lower range. This difference may be caused by the heat loss to the support fiber of 
the fuel droplet. It can be noted that the Kumagai et al. conducted experiments 
using unsupported droplet and which may not have heat losses to the support fiber 
as compared to the data of Xu et al. [79]. Hence it is expected that the present 
prediction can match with the experimental data of Kumagai. This small difference 
may be due to consideration of gray gas radiation model and the finite rate 
chemistry. In contrast, the present predicted data at high temperature does not 
match well with the experimental data of Xu et al. [79]. However its trend is 
qualitatively similar to the experimental data. 



Fig.5.17: Transient variation of flame temperature for (a) ambient room 
temperature and (b) ambient temperature = 1093 K. 

In order to investigate the present anomalies between the predicted data and 
the experimental data of Xu et al. [79] at the high temperature, the variation of the 
flame temperature during droplet burning is plotted in Fig. 5.17 (b) with different 
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initial droplet diameter. It can be observed from the Fig. 5.17(b) that at high 
ambient temperature the flame temperature become very high almost from 2800 K 
to 3600 K for droplet diameter of 1.0 to 2.0 mm. The predicted flame temperature 
seems to be quite high which may not be possible in the actual situation due to 
occurrence of endothermic dissociation reaction. However the present model does 
not contain the dissociation reactions which need to be included to predict this 
aspect. 

The variations of flame temperature with time for three different initial 
droplet diameters are plotted in Fig 5.17 (a) at ambient temperature. It can be 
observed from Fig. 5.13 (a) that, the flame temperature is found to be decreased 
with the droplet diameter. In contrast the flame temperature increases with the 
diameter at high ambient temperature of 1093 K as shown in Fig. 5.17(b). This 
increase in flame temperature is expected to enhance the burning rate of the fuel 
droplet, which was discussed earlier at length. 

5.2. Multicomponent Fuel Droplet Combustion 

In modem days, practical fuels, used in combustion devices are blends of several 
hydrocarbons with different chemical and physical properties. The properties of 
such kind of fuels can be approximated with two hydrocarbons with different 
volatilities. Hence, the present study is carried out with a binary fuel consists of n- 
Heptane and n-Decane. Special efforts are made on the transient analysis of 
temperature, species mass fraction profile for both the gas phase and liquid phase 
by varying the initial fuel compositions. Transient variations of the surface 
parameters are very important to study in the case of multicomponent fuel droplet 
combustion. It has been considered in the literature that internal circulation takes 
place inside the droplet. Hence, the distillation curve model is being used in the 
present study by assuming the occurrence of proper mixing of fuels inside the 
droplet. The following subsections are devoted to discuss the results of 
multicomponent fuel droplet combustion. 


dVd 


66 


5.2.1. Effect of Composition on Burning Parameters 

In the present section the influences of the compositions of the binary fuel on the 
droplet burning are discussed. The droplet diameter square is plotted in Fig. 5.18 (a) 
with the time with variation in fuel compositions. It can be observed from this 
figure that burning rate and total time required for binary fuel droplet combustion 
vary as the composition of fuel components are changed. Average burning rate 
constant is found to be enhanced with the increase in high volatile fuel composition 
as per expectation. Although the pure fuel droplet combustion exhibit classical cf 
law behavior but results of binary fuel droplet with variation in compositions does 
not seem to follow the cf law behavior. 



t/d] (sec/mm^) t/d^ (sec/mmT 


Fig.5.18: Effect of initial compositions of n-heptane and n-decane on (a) law and 
(b) Burning rate constant Kb- 

To understand the burning rate behavior in details, the burning rate constant, 
Kb with the variation in fuel compositions is plotted with the time in Fig. 5.18 (b). A 
typical behavior corresponding to the binary fuel with compositions of 50% n- 
Heptane-50% n-Decane can be observed from this figure. Initially the burning rate 
constant, Kb increases very sharply in the initial preheat period. During this period 
droplet surface ternperature also increases and attains a stable value. The 
vaporization of the high volatile substances is likely to be more than the less 
volatile fuel. It can be observed from the Fig. 5.18 (b) that the burning rate attains 
almost a peak value after the initial preheat period during which the vaporization of 
high volatile fuel is likely to occur at a faster rate. After certain time the high 
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volatile fuel will be depleted quickly. Then the burning rate constant starts 
decreasing and attains its minimum value. During this time both the components of 
the binary fuel droplet attains the same gasification rate as shown in Fig. 5.19. It 
can be noted that the surface temperature attains an inflection point at this time. 
Subsequently, the bummg rate starts increasing due to the higher gasification rate of 
n-Decane fuel as shown in Fig. 5.19. It can also be observed in Fig. 5.19 that the 
surface temperature increases monotonically towards the boiling point temperature. 



Fig.5.19: Transient variations of the Surface Parameters for n-Heptane and n- 
Decane fuel mixture of compositions ratio 50%-50%. 

In case of briny fuel droplet with composition of 90% n-Heptane- 10% n- 
Decane, the burning rate constant attains its minimum value towards the end of the 
burning as shown in Fig. 5.20 (a). As concentration of the high volatile n-Heptane 
fuel is higher than the low volatile n-Decane fuel, it is expected that the burning of 
n-Heptane fuel will dominant for longer period of time maintaining a lower 
temperature at its surface. As a result, n-Decane fuel can play an important role 
during the combustion of binary fuel after the complete consumption of n-Heptane 
fuel, which can be seen in Fig. 5. 20 (a). Thus the surface temperature increases at a 
faster rate' towards the boiling point temperature of n-Decane fuel. As a result the 
burning rate drops down as shown in Fig. 5.20 (a). 

However, in case of binary fuel, composed of 10% n-Heptane-90% decane 
n-Heptane gets vaporizes initially and gets depleted within short period of time. The 
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gasification rate of the heptane is found to be decreased with time as the percentage 
of n-Heptane is decreased in binary fuel. However the vaporization rate of n- 
Decane is predicted to be higher than the n-Heptane due to higher percentage of n- 
Decane in binary fuel. The temperature at droplet surface keeps increasing as shown 
in Fig. 5.20 (b) till n-Heptane fuel gets vaporized completely and then it maintains 
almost a constant value. As a result the burning rate increases with the time. 
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Fig.5.20: Transient variations of the Surface Parameters for n-Heptane and n- 
Decane fuel mixture of compositions ratio of (a) 90%-10% and (b) 10%-90%, 



Fig.5.21: Effect of initial compositions of n-heptane and n-decane on (a) droplet 
surface temperature and (b) variation in the droplet surface temperature with time. 
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The transient variations of droplet surface temperature are plotted in Fig. 
5.21 (a). Depending on its compositions, the droplet surface temperature increases 
from a stable value towards the higher boiling point temperature (n-Decane). The 
temperature variation rate at the droplet surface is plotted in Fig. 5.21(b) with 
normalized time for binary fuel with representative composition. It can be observed 
that the temperature variation rate attains a maximum value after certain time from 
the initiation of droplet burning. Its peak moves forward and is increased in 
magnitude with the increase in the high volatility component in the binary fuel 
droplet. Interestingly, the temperature variation rate becomes steeper with increase 
in the high volatile component indicating that surface temperature increases at a 
faster rate towards the end of the burning. 

The effects of the fuel compositions on the droplet center temperature are 
plotted in Fig.5.22 (a) with normalized time. For pure n-Heptane and n-Decane the 
droplet center temperature increases monotonically towards their respective boiling 
point temperature. The temperature variation rate at the center of droplet is plotted 
in Fig. 5.22 (b) for representative fuel compositions. It can be noted that the slope 
of temperature profile for the binary fuel droplet gets varied due to change in its 
surface temperature. A sharp change can be noticed towards end of burning for fuel 
composition of 90% n-Heptane- 10% n-Decane. However for the fuel composition 
of 50% n-Heptane-50?4 n-Decane, the change in the slope is occurred almost in the 
middle of the droplet burning. 



Fig.5.22: Effect of initial compositions of n-heptane and n-decane on (a) droplet 
center temperature and (b) variation in the droplet center temperature with time. 
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Furthermore, the fractional mass gasification rate for n-Heptane and n- 
Decane are plotted in Fig. 5.23 (a) and (b) respectively. The fractional mass 
gasification variations rates for n-Heptane and n-Decane are plotted in Fig. 5.24 (a) 
and 5.24(b) respectively. It can be observed from the Fig. 5.23(a) that the fractional 
mass gasification rate for n-Heptane decreases with its concentration. In contrast, 
the fractional mass gasification rate for n-Decane increases with its concentration. 
Similar features can be also observed in the temporal variations of mass fractions as 
plotted in Fig. 5.25 (a) and (b). 




Fig.5.23: Effect of initial compositions of n-heptane and n-decane on fractional 
mass gasification rate for (a) n-Heptane and (b) n-Decane. 



Fig.5.24: Effect of initial compositions of n-heptane and n-decane on rate of change 
of fractional mass gasification rate for (a) n-Heptane and (b) n-Decane. 
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Since the boiling point temperature of n-Heptane is less than n-Decane, the 
gasification of n-Heptane is expected to occur at higher temperature. Initially more 
amount of n-Heptane is likely to be gasified and its mass fraction inside the droplet 
gets decreased at faster rate. As a result, the gasification rate is decreased further 
leading to the complete depletion of n-Heptane fuel. Subsequently the n-Decane 
fuel gets gasified at higher surface temperature of the droplet. 
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Fig.5.25: Effect of initial compositions of n-heptane and n-decane on liquid phase 
mass fractional profiles for (a) n-Heptane and (b) n-Decane. 
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Chapter 6 
CONCLUSIONS 


6.1. Conclusions 

In this work the spherical symmetric combustion of single isolated fuel droplet has 
been investigated both experimentally and numerically. Experimental investigation 
of droplet combustion for both single and binary fuel has also been performed under 
normal gravity condition at ambient pressure and temperature. In order to simulate 
isolated droplet combustion, one dimensional, unsteady governing equations are 
solved numerically using Finite Difference Method. In present computation a single 
step chemistry model is employed. In this CFD code variable properties and gray 
gas model for radiation heat transfer have been considered. The effect of ambient 
oxygen level and effect of radiation heat transfer on the burning parameter at 
different temperatures have been studied extensively. The influence of initial 
droplet diameter on the droplet burning rate is also investigated. Besides this, the 
multicomponent fuel droplet combustion is simulated using Distillation Curve 
Model. The major conclusions of the present study are as follows; 

Conclusions from the Experimental Data: 

1. Both the single and multicomponent fuel droplet combustion follow the 
classical </law behavior under normal gravity at ambient pressure and 
temperature. 

2. Initially flame envelops increases and subsequently diminishes with time. 

3. The slope of the cf law curve known as the burning rate constant, Kb are 
found to be almost constant for both single and multicomponent fuel droplet 
combustion. 

4. The burning rate constant, Kb for both n-Heptane and n-Decane is found to 
be 1.12 and 1.03 respectively. For binary fuel droplet with composition of 
50% n-Heptane and 50% n-Decane, the burning rate constant, Kb is found to 


be 1.1. 
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Conclusions from the Numerical Studies: 

1 • Th® burning rate of n-Heptane fuel droplet at normal ambient conditions is 
predicted to be higher than the existing experimental data. This 
discrepancies is attributed the considerable radiation heat transfer from the 
high temperature flame zone to its surrounding. Hence a suitable radiation 
model is thought to be included in the present code. 

2. Initially, the flame standoff distance is found to be increased and 
subsequently diminishes with time. This phenomenon is commensurate to 
our experimental data. The flame stand-off ratio varies vdthin 11-28. But in 
the literature the experimental flame standoff ratio is reported to be within 
4-10. To predict the flame standoff ratio accurately radiation heat transfer 
model and multi step kinetics chemistry model may be required. 

3. The present numerical prediction indicates that the burning rate gets 
enhanced with increase in ambient oxygen concentration level. The increase 
in the ambient oxygen concentration level is found to enhance the reaction 
rate, resulting in an increase in flame temperature. The heat transfer from 
the flame zone to the droplet surface causes the enhancement of gasification 
rate. 

4. The flame stand-off ratio is found to be increased with the ambient oxygen 
level. This is in contradiction with the classical quasi-steady theory. These 
anomalies may be due to fact that the classical theory does not take account 
of the gas phase velocity, which is the determining factor in predicting the 
flame standoff ratio. 

5. The gray gas model for radiation heat transfer is incorporated in the present 
code. As a result, the numerical prediction improves considerably in 
accordance with the experimental data. The predicted flame temperature is 
found to be closer to the experimental data with the inclusion of radiation 
heat transfer model. The burning rate data also matches well with the 
existing experimental data. 



It is found out that the radiation heat transfer plays a predominant role in 
determining the burning rate particularly for bigger droplet. As the droplet 
diameter decreases the heat loss from the flame surface decreases leading to 
the rise in burning rate. 

By varying the initial droplet diameter, it is found out that burning rate 
diminishes at ambient temperature. In contrast, an opposite trend in burning 
rate is observed at high ambient temperature. The present predicted data at 
low ambient temperature matches well with the experimental data of Xu et 
al [79]. However, at high ambient temperature the present prediction seems 
to match only qualitatively with the experimental data. The prediction may 
be improved by inclusion of dissociation reactions in the present model. 

During the burning of multicomponent fuel droplet, a slight variation from 
the classical cf law is observed. The burning rate and total time for burning 
varies with the composition of the droplet. 

Depending on the composition of fuel, the variation in the burning rate 
constant is observed. During the course of burning it has been observed that 
the droplet surface temperature is increased from a low value to a high value 
corresponding to the boiling point temperature of respective constituents m 
the binary fuel. 

Initially the high volatile fuel is found to be getting consumed at a faster rate 
in comparison to low volatile fuel. It has been observed that the fractional 
mass gasification rate for n-Heptane decreases with its concentration. In 
contrast, the fractional mass gasification rate for n-Decane is found to be 
increased with its concentration. 
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6.2. Recommendations for Further Studies 

During this work it has been observed that there are many issues that need further 
investigation. Some of them are listed below. 

i) Multi-step chemistry model is to be incorporated in the computer code to 
improve its prediction. 

ii) Droplet combustion at super critical pressure. 

iii) The interaction between the droplets can be investigated. 

iv) The model for multicomponent fuel droplet combustion can be improved 
further. 

v) Soot production and soot radiation study may be carried out. 

vi) High temperature dissociation effect can be considered. 
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Appendix A 

Liquid and Gas Phase Thermal and Transport Properties 

Calculation 

1. Gas Phase Properties 


(a) Density: 

The density of the gas mixture, p, has been obtained from the equation of state: 


P = 


RJ 


(A.1) 


where, P is pressure; is mass fraction of species i ; m. is molecular weight of 
species i , and is universal gas constant. 


(b) Specific Heat 

The specific heat of the gas mixture has been estimated assuming it to be an ideal 
gas mixture of CO 2 , H 2 O, O 2 , N 2 and fuel vapor. The specific heat, Cp, of mixture 
can be calculated by using the following relation. 

(A.2) 

/=1 

The specific heat of component can be calculated from the curve-fit 
correlation. The data is taken from Sandia database and Toulokian et al. [72]. Then 
verify the data with Virgoftik [74]. 

^=(c.+c,r,+c//+c,r/+c.j;‘) for cch.ftO.Oj.N, (a.3) 
Cp^ (J lkgK) = (Q + C.T; + C// + CjT/ + c//)x 41 87 for fuel vapors (A.4) 

The value of the coefficient Cq , C, , C 2 , C 3 and C 4 are listed in table A. 1 . The 
plot of the Specific Heat with time is shown in Fig. A.l. 
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Table A.1: Constants for the calculation of constant pressure specific heat. 


Species 

Temp. Range(K) 

Co 

c, 

C2 

C3 

C4 

CvH.e 

300 - 800 

0.0254671 

1.37666e-3 

-4.0004e-7 

-1.1331e-10 

0.0 

C 7 H ,6 

800-2000 

0.0982018 

1.30398e-3 

-5.92804e-7 

1.0105e-10 

0.0 

C,oH22 

300 - 780 

0.1273 

8.9679e-4 

4.2932e-7 

-5-1068e-10 

0,0 

C,oH22 

780 - 2000 

0.0188191 

1.5129e-3 

-6.7117e-7 

1.1236e-10 

0.0 

O2 



1.483e-3 

-7.579e-7 

2.0947e-10 

-2.167e-14 

O2 



-2.9967e-3 

9.847e-6 

-9.6813e-9 

3.2437e-12 

0 

u 

200 - 1000 

3.8574 

4.414e-3 

-2.2148e-6 

5.2349e-10 

-4.7208e-14 

0 

0 

1000-3500 

2.3567 

8.9846e-3 

-7.1235e-6 

2.459e-9 

-1.437e-13 

H20 

200 - 1000 

3.03399 

2.1769e-3 

-1.640e-7 

-9.704e-l 1 

1.682e-14 

H20 

1000 - 3500 

4.1986 

-2.0364e-3 

6.5204e-6 

-5.4879e-9 

1.772e-12 

N2 

300 - 1000 

2.9266 

1.4879e-3 

-5.684e-7 

L0097e-10 

-6.7533e-15 

N2 

1000-5000 

3.2986 

1.40824e-3 

-3.9632e-6 

5.6415e-9 

-2.4448e-12 



Fig.A.l: The variation of Specific Heat of Different Species with Gas Phase 
Temperature. 
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(c) Thermal Conductivity 


Gas mixture thermal conductivity has been calculated by using Wassiljewa equation 
(1904) and Masomn and Sexana modification (1958), which has been 
recommended for both polar and non-polar gas mixtures (Reid et al. [59]). 

(A.5) 


/=:1 




j=i 

where, is the thermal conductivity of the gas mixture, thermal conductivity of 
pure componet i . Xj is the mole fraction of the component i . Ay can be found out 
from the Mason and Saxena (1958) modification. 



(A.6) 


where, m, , trij are the molecular weight of components i and j . 


Table A.2: Constants for the calculation of Thermal conductivity (W/m s). 


Species 

Q 

C, 


Q 

Q 1 

C 7 H 16 

4.7317e-3 

-4.3918e-8 

2.8322e-7 

-1.8841e-10 

3.8818e-14 

C 10 H 22 

3.9817e-3 

-3.0918e-8 

2.8352e-7 

-1.8841e-10 

3.882e-14 

O 2 

1.3708e-3 

9.2262e-5 

! 

-2.8054e-8 

6.3322e-12 

-2.299e-16 

CO 2 

1.2528e-3 

1.8145e-5" 

1.4728e-7 

-1.3644e-10 

3.7908e-14 

H 2 O 

2.0878e-2 

-9.4108e-5 

3.6958e-7 

-2.7724e-10 

6.7850e-14 

N 2 

3.7894e-3 

8.1514e-5 

-2.339e-8 

-3.2294e-12 

4.3092e-15 


The thermal conductivity of the pure component can be obtained by taking 
data from Toulokian et al. [72] and finding out cxirve-fit correlation. 


= (Q + c,r, +C 3 r/ +C,7’/) (A.7) 

The value of the coefficient Q, C, ,^,€ 3 and C, for thermal conductivity 
are listed in table A.2. The variations of thermal conductivity with the temperature 
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for different species have been plotted in Fig. A.2. 



Fig. A. 2: The variation of Thermal Conductivity of Different Species with Gas 
Phase Temperature. 


(d) Diffusion Coefficient 

Diffusion coefficient of any component, £>, , in the mixture can be calculated as 
follows (Wilke, 1950) 


A = 


1 -^, 

A, A, 


(A.8) 


The binary gas diffusion, D^, has been calculated from the following 
expression ( Hirschfelder [22]) 


A, =262.8x10' 




(m, + ) 

2ffj, 


PayCl/y 


(A.9) 


where, P is the pressure in bar, and /Mj are the molecular weight of the 
components. cTyis the collision diameter in A, cr^ =0.5(cr; +<7^). Q^^is the 
collision integral, defined as Perry [56] 
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=(44.54x7;‘-^’'>’ +1.911x7:;-'^”)“' (A.IO) 

41 I « JT 

where, T^j = ^jT^ Tj , aiid 7]. = —.The value of and Ik taken from Sandia 

Go / /C 

Database is tabulated in Table A.3 . 


Table A.3: Constants for the calculation of Transportation Data 


Species 

C7H,6 

C 10 H 22 

O 2 

CO 2 

H 2 O 

N 2 

cr, (A) 

6.253 

7.229 

3.458 

3.763 

2.605 

3.621 

Go I k (K) 

459.6 

600.6 

107.4 

244 

572.4 

97.530 


2. Liquid Phase Properties 

The Liquid phase properties, density, specific heat and thermal conductivity have 
been calculated using curve-fit correlations given in Touloukian et al. [72] and 
checked from the experimental data given in Vargafitik, [74], 

(a) Density 

For calculating the density for n-Heptane and n-Decane following correlations are 
used respectively. 

= 720.0-0.892(7; -253.0) forn-Heptane (A.11) 

/ 9 ;(itg/m'') = 923.4-0.5587xr-0.0003xr^ forn-Decane (A.12) 

(b) Specific Heat 

Correlation used for specific heat of liquid is 

C,j {J ! kg K) = (Co + C,7; + ) x 41 87 (A. 1 3) 

The values of the coefficients are as follows: 

For n-Heptane this correlation is valid between 230K to 480K and 
Co =0.454958, C, =-3.46073 xl0'\ C^ =2. 1482x1 0"^ Cj =-2.48653x10'^ 
For n-Decane this correlation is valid between 283 K to 493 K and 
C, =0.581718, C,= -8.44526x10-*, C, =1.21946x10^. C, =3.37247x10-’ 
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(c) Thermal Conductivity 

The following correlations are used to obtain the thermal conductivity for n- 
Heptane and n-Decane respectively 

X,{JlmK) = (466.94 - 0.54892 x 7) ) x 4. 1 87 x 1 O'" (A. 1 4) 

;i, (J / m is:) = (564.599 - 0.78 1 586 X 2; ) X 4. 1 87 X 1 0^ (A. 1 5) 

These correlations are gives substantially accurate values from 290 to 370 K and 
230 K to 450 K respectively, and outside these range uncertainty increases. As our 
domain of interests is within these, these correlations are acceptable. 

(d) Latent Heat of Vaporization 

Latent heat of evaporation is a function of the surface temperature, , at low 

pressure, is approximated by the following equation taken from Perry [55] and 
checked with V argaftik [74] . 

UJIkg)=C,(Tc-T,f (A.16) 

where, 7^ is the critical temperature, for h-Heptane and n-Decane 540.2K and 
617.7K respectively. Cf) and C,are the two coefficients. For n- 
HeptaneCg =55.53 xl0\ C, =0.3436 and for n-Decane Q =36.02x10^, 
C, =0.39797. 

3. Radiative Properties (Planck Mean Absorption Coefficient) 

The Planck Mean Absorption Coefficient of CO 2 and H 2 O, and , have been 

determined by the curve-fit data taken from RADCAL program 1993 [19] as sited 
in http://www.ca.sandia.gov/TNF/radiation.html. These curve fits were generated 
for temperatures between 300K and 2500K and may be very inaccurate outside this 
range. 

= Co+c,x(i000/7;)+c,x(i000/rp^+C3x(i000/r,)^ 

+C4 X (1 000 / Tg )“ + Q X (1 000 / Tg )' 

The values of the Coefficients are given in Table A.4. The Planck Mean Absorption 
Coefficients for CO2 and H 2 O are plotted in Fig. A.4. 
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'Table A.4: Constants for the Calculation of 
Coefficient 


Planck Mean Absorption 


Coefficients 

Co 

Cl 

C 2 

Ci 

C 4 

Cs 

CO 2 

18.741 

-121.310 

273.5 

-194.05 

56.310 

-5.8169 

H 2 O 

-0.23093 

-1.12390 

9.41530 

-2.99880 

0.51382 ■ 

-1.8684e-5 



Fig.A.3: Planck Mean Absorption Coefficients for H 2 O and CO 2 . 
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